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SECTION  I 


INTRODUCTION 

Flows  in  which  a  separated  turbulent  shear  layer  reat¬ 
taches  after  formation  exist  in  many  problems  of  engineering 
importance  including  external  flows  over  structures  and 
internal  flows  with  sudden  expansions.  Axisymmetr ic  sudden 
expansions  are  of  particular  interest  because  the  region  of 
recirculation  can  be  used  to  advantage  as  a  flame  holder  for 
a  dump  combustor .  In  spite  of  this  importance  and  the  large 
body  of  experimental  data  available  for  such  flows,  there  is 
still  an  incomplete  understanding  of  the  problem.  Part  of 
the  difficulty  lies  in  the  fact  that  these  flows  are  highly 
turbulent  with  frequent  velocity  reversals  which  severely 
limits  the  quality  of  measurements  made  with  conventional 
techniques  such  as  pitot  tube  and  hot  wire  anemometry. 
Reacting  flows  complicate  the  measurement  problems  even 
further.  With  the  laser  Doppler  velocimeter  (LDV) ,  however, 
valuable  information  about  such  flows  can  be  gained.  This 
instrument  has  many  desirable  qualities:  no  physical  probe 
to  intrude  on  the  flow,  high  spatial  and  temporal  resolu¬ 
tion,  linear  response,  and  the  capability  of  determining  the 
direction  of  the  velocity  being  measured.  A  complete 
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description  of  the  principles  of  laser  Doppler  velocimeters 
is  given  by  Stevenson  [1).  * 

Although  the  LDV  is  ideally  suited  for  flow  measure¬ 
ments  m  highly  turbulent  mixing  flows  with  flow  reversal, 
it  is  Known  that  certain  bias  errors  can  occur  [2].  There¬ 
fore  one  objective  of  the  present  study  was  to  further  ver¬ 
ify  a  previously  developed  experimental  technique  for  elim¬ 
inating  velocity  bias  [3,4,5].  In  the  present  investigation 
three  complete  sets  of  LDV  measurements  were  made  in  an 
axisymmetric  sudden  expansion.  The  three  data  sets  con¬ 
sisted  of  mean  streamwise  velocity  and  turbulence  intensity 
measurements,  two  of  which  were  made  in  isothermal  (cold) 
flow  (unbiased  and  biased)  and  one  of  which  was  made  in  a 
reacting  (hot)  flow  (biased).  (Unbiased  hot  flow  measure¬ 
ments  could  not  be  made  because  of  seeder  limitations.)  The 
biased  and  unbiased  cold  flow  measurements  were  compared  to 
determine  the  effect  of  velocity  bias  on  the  flow  measure¬ 
ments.  The  hot  and  cold  flow  measurements  were  compared  to 
determine  the  effect  of  combustion  on  the  turbulent  struc¬ 
ture  of  the  flow  field. 

A  review  of  the  general  flow  characteristics  of  axisym¬ 
metric  sudden  expansions  together  with  a  review  of  LDV  meas¬ 
urements  in  reacting  flow  fields  are  presented  in  Section 
II.  A  discussion  of  measurement  errors  which  occur  in  laser 

1 .  UumEer s  in  brackets  refer  to  references  listed  at  the 
end  of  the  report. 
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\  velocimetiy  is  also  given  in  Section  II.  Sections  (II,  IV 

and  V  describe  the  apparatus,  techniques  and  results  for  the 
experimental  measurements.  A  description  of  the  numerical 
code  and  a  comparison  of  numerical  predictions  to  the  isoth¬ 
ermal  experimental  results  is  given  in  Section  VI.  Conclu¬ 
sions  and  recommendations  are  presented  in  Section  VII. 
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SECTION  II 


LITERATURE  REVIEW 


1.  INTRODUCTION 

The  general  flow  characteristics  of  axisymmetr ic  sudden 
expansions  are  of  great  engineering  interest  and  have  been 
the  focus  of  many  recent  experimental  and  analytical  inves¬ 
tigations.  One  purpose  of  this  section  is  to  describe  the 
flow  field  of  an  axisymmetr ic  sudden  expansion  and  identify 
some  of  the  relevant  mechanisms  responsible  for  the  charac¬ 
ter  of  the  flow  field.  A  review  of  the  recent  literature 
pertaining  to  axisymmetric  sudden  expansions  is  also 
included. 

Combustion  significantly  changes  the  character  of  the 
flowfield  by  introducing  steep  temperature  gradients.  The 
effects  of  combustion  on  the  flow  field  and,  m  particular, 
on  the  turbulence  in  the  flow  has  also  been  of  great 
interest.  Much  work  has  been  conducted  in  an  effort  to 
understand  the  turbulence-combustion  interaction.  Some 
recent  literature  dealing  with  LDV  measurements  in  reacting 
flows  will  also  be  reviewed.  Finally,  a  discussion  of 
errors  associated  with  LDV  measurements  will  be  presented. 
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2.  AX I SYMMETRIC  SUDDEN  EXPANSIONS 

a.  PAST  EXPERIMENTS 

Several  studies,  both  experimental  and  computational, 
have  been  conducted  to  determine  the  effect  of  such  parame¬ 
ters  as  inlet  flow  characteristics  on  the  shear  layer  in 
axisymmetric  sudden  expansions.  Experimenters  have  used  a 
variety  of  techniques  to  study  the  velocity  field  in  flows 
of  this  type  including  flow  visualization,  hot  wire  and  hot 
film  anemometry,  and,  most  recently,  laser  Doppler  velo- 
cimetry.  Air  and  water  have  been  the  basic  fluids  used, 
with  only  limited  work  conducted  with  reacting  gases.  Rey¬ 
nolds  numbers  based  on  average  inlet  velocity  and  inlet 

7 

diameter  typically  range  from  10  to  10  .  The  boundary  layer 
at  the  point  of  separation  has  been  either  laminar  or  tur¬ 
bulent. 

Table  1  is  a  compilation  of  recent  investigations  in 
axisymmetric  sudden  expansions.  This  is  certainly  not  a 
complete  list,  but  includes  those  studies  pertinent  to  the 
discussion.  Recent  reviews  of  the  related  problem  of  two- 
dimensional  step  flows  are  given  by  Brenner,  et  al.  (3]  and 
by  Eaton  and  Johnston  [6,7]. 
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TABLE  1.  DETAILS  OF  PREVIOUS  AXISYMMETRIC  ANNULAR  STEP  INVESTIGATIONS. 


b.  SEPARATED  PLOW  REGIONS 


The  separated  flow  downstream  of  a  sudden  expansion  can 
be  divided  into  two  regions:  1)  the  mixing  layer  region  and 
2)  the  relaxation  region  as  shown  in  Pigure  1.  The  mixing 
layer  region  includes  the  flow  from  separation  to  reattach¬ 
ment.  The  mixing  region  may  be  further  subdivided  into  four 
flow  regions,  namely  I)  The  secondary  recirculation  zone 
(characterized  by  one  or  more  vortices  thought  to  be  rotat¬ 
ing  about  an  axis  perpendicular  to  the  radial  axis);  II) 
The  primary  recirculation  zone  consisting  of  trapped  eddies 
which  rotate  clockwise  (in  the  bottom  recirculation  zone)  in 
a  meridional  plane;  III)  The  reattachment  region  where 
bifurcation  of  the  shear  layer  occurs  and  part  of  the  flow 
is  deflected  upstream  into  Region  1 1  to  supply  entrainment; 
and  IV)  A  curved  free  shear  layer  with  a  varying  velocity 
deficit  across  it,  characterized  by  high  turbulence  and 
large  intermittent  eddies  that  promote  the  mixing  process. 

The  relaxation  region  begins  at  flow  reattachment  and 
ends  where  full  recovery  of  the  turbulent  boundary  layer 
occurs.  The  low  turbulence  core  and  high  turbulence  shear 
layer  lose  their  identity  in  this  region  and  the  mean  velo¬ 
city  profile  approaches  that  expected  for  a  fully  turbulent 
pipe  flow.  Studies  have  defined  the  turbulent  structure  of 
Regions  II  and  III,  but  the  structure  in  Region  I  has  yet  to 
be  conclusively  identified  in  axisymmetr ic  sudden 
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expansions.  This  is  probably  due  to  the  low  velocities, 
typically  1%  of  the  inlet  velocity,  present  in  this  region. 

One  would  expect  symmetric  flow  patterns  in  axisym- 
metric  sudden  expansions  for  the  following  reason.  An  asym¬ 
metric  flow  would  create  an  asymmetric  pressure  distribution 
in  the  separation  region.  This  would  be  offset  by  a  redis¬ 
tribution  of  pressure  within  this  region  leading  to  a  sym¬ 
metric  flowfield.  In  a  laminar'*'  flow  study.  Macagno  and 
Hung  [8]  showed  that  symmetric  flow  patterns  are  maintained 
over  a  wide  range  of  Reynolds  numbers  for  axisymmetr ic  sud¬ 
den  expansions.  Zemmanic  and  Dougall  [9]  did  observe  an 

asymmetric  flow  pattern  (determined  by  heat  transfer  raeas- 

2 

urements)  for  turbulent  flow  of  air  in  an  axisymmetr ic  sud¬ 
den  expansion.  The  extent  of  asymmetry  was  small  compared 
with  that  which  occurs  in  a  plane  expansion,  however.  This 
is  the  only  reference  found  which  noted  any  asymmetry  in  an 
axisymmetr ic  sudden  expansion  flow. 


c .  RE ATTACHMENT  LENGTH 

The  streamwise  distance  between  the  point  of  separation 
and  reattachment  for  axisymmetr ic  sudden  expansions  has  been 
the  subject  of  many  investigations.  Functional  relation¬ 
ships  between  reattachment  length  and  Reynolds  number,  step 
1-  i-e.  a  flow  with  a  laminar  boundary  layer  at  separation 

2.  i.e.  a  flow  with  a  turbulent  boundary  layer  at 
separation 
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height,  H  (or  area  ratio,  AR) ,  and  inlet  flow  conditions 
have  been  postulated. 

Back  and  Roscke  [10]  performed  dye  studies  in  an 
axisymmetr ic  sudden  expansion  (oriented  horizontally)  using 
water  as  the  working  fluid  for  a  Reynolds  number  (based  on 
inlet  pipe  diameter  and  average  inlet  velocity)  ranging  from 
20  to  4200.  The  purpose  of  the  investigation  was  to  study 
the  effect  of  inlet  Reynolds  number  on  shear  layer  growth 
and  t eat.tachment  length.  Their  results  showed  that  laminar , 

3 

transitional,  and  turbulent  separated  flow  occurred  within 
the  Reynolds  number  range  of  their  study  with  reattachmemt 
length  varying  significantly.  Reattachment  length  had  a 
maximum  value  of  25  step  heights  for  an  inlet  Reynolds 
number  of  290.  For  Reynolds  numbers  greater  than  290  lam¬ 
inar  instabilities  became  visible  in  the  inlet  boundary 
layer  and  the  reattachment  length  decreased  very  rapidly  to 
about  seven  step  heights  at  Re: 1000.  The  reattachment  length 
then  slowly  increased  to  a  nearly  constant  value  of  approxi¬ 
mately  nine  step  heights  for  Reynolds  numbers  greater  than 
3000. 


Freeman  [11]  measured  streamwise  velocities  of  water  in 

the  turbulent  flow  field  of  an  ax i symmetric  sudden  expansion 

(oriented  vertically)  using  an  LDV  system.  A  frequency 

locked  tracking  processor  along  with  a  digital  voltmeter  and 

~~T~.  i .  e .  a  flow  with  a  laminar  or  transit  loning  boundary 
layer  at  separation  which  becomes  turbulent  by 
reattachment. 
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rms  meter  were  used  to  measure  mean  and  fluctuating  veloci 


ties,  respectively.  Based  on  experimental  stream  function 
contours,  reattachment  was  found  to  occur  at  8.7  step 
heights.  Maximum  negative  recirculation  velocities  were 
aproximately  10%  of  the  inlet  centerline  velocity. 

Moon  and  Rudinger  [12],  m  a  turbulent  axisymmetric 
sudden  expansion  experiment  using  air,  found  the  reattach¬ 
ment  length  to  be  between  8  and  9  step  heights.  They 
located  the  reattachment  point  by  interpolating  mean  velo¬ 
city  profiles  obtained  with  an  LDV  at  several  cross- 
sectional  planes.  They  found  that  the  flow  field  downstream 
of  the  sudden  expansion  was  symmetric,  thus  contradicting 
2emmanic  and  Dougall  [9}.  They  also  concluded  that  reat¬ 
tachment  length  had  no  functional  dependence  on  Reynolds 
number  for  turbulent  flows. 

An  experimental  study  was  performed  by  Drewry  [13]  on 
the  axisymmetric  sudden  expansion  in  ramjet  combustor  models 
(cold  flow)  using  flow  visualization  techniques,  wall  static 
pressure  measurements,  and  gas  sampling.  Drewry  found  that 
reattachment  length  varied  linearly  with  step  height.  The 
reattachment.  point  occurred  between  7.9  and  9.2  step  heights 
for  all  geometries  tested.  He  also  found  the  flow  to  be 
symmetric  with  slight  circumferential  flow  (swirl)  down¬ 
stream  of  the  step.  No  explanation  was  given  for  this  swirl 
component,  but  it  may  have  been  introduced  by  the  inlet  air 
facility. 
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Kankovi  and  Page  [14],  using  static  pressure  taps  and 
hot  wire  anemometry,  found  that  reattachment  occurred  at 
about  8  step  heights  downstream  of  the  sudden  expansion. 
They  also  identified  the  presence  of  a  weak  secondary  recir¬ 
culation  zone  (Region  I)  located  within  one  step  height  of 
the  sudden  expansion. 

Stevenson,  et  al.  [15],  in  a  turbulent  axisymmetric 
sudden  expansion  experiment  with  air,  measured  streamwise 
and  tangential  mean  velocities  and  turbulence  intensities 
using  a  one  component  LDV  system.  They  also  derived  Rey¬ 
nolds  stress  correlations  from  independent  measurements  of 
velocity  components.  Logan's  method  [16J  was  used  to  derive 
the  Reynolds  stress  from  the  individual  component  measure¬ 
ments.  They  found  that  i ^attachment  occurred  at  7.9  step 
heights  downstream  of  the  sudden  expansion.  The  flow  field 
was  also  found  to  be  symmetric. 

A  review  of  the  studies  described  above  show  that  the 
observed  reattachment  length  for  axisymmetric  sudden  expan¬ 
sions  varies  between  7  and  10  step  heights  for  turbulent 
inlet  flow  conditions.  Although  this  was  noted  by  the 
experimenters,  little  has  been  done  to  determine  the  reason 
for  the  variation.  Eaton,  et  al.  [17]  suggested  that  part 
of  the  variation  of  reattachment  length  for  similar 
geometries  and  inlet  flow  conditions  may  be  the  result  of 
the  recirculation  zone  slowly  growing  and  shrinking  causing 
the  reattachment  streamline  to  oscillate.  The  entrainment 
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cate  balances  the  backflow  rate,  but  only  in  the  mean,  not 
instantaneously.  Limited  observations  suggest  the  oscilla¬ 
tions  to  be  of  relatively  low  frequency  (on  the  order  of  50 
Hz).  Although  there  is  some  experimental  evidence  to  sup¬ 
port  the  unsteady  reattachment  postulate,  the  phenomena 
needs  to  be  investigated  in  much  more  detail.  Keuhn  [16] 
suggested  that  part  of  the  reattachment  length  variation  can 
be  attributed  to  differences  m  inlet  conditions,  but  that 
most  may  be  due  to  an  adverse  pressure  gradient  effect.  He 
showed  that  superimposing  different  pressure  gradients  on 
the  sudden  expansion  flow  field,  while  keeping  inlet  condi¬ 
tions  constant,  led  to  large  variations  in  reattachment 
length.  The  reattachment  length  for  laminar  and  transi¬ 
tional  flows  has  a  strong  dependence  on  the  state  of  the 
boundary  layer  and  the  turbulence  level  in  the  free  stream 
[10].  For  turbulent  flows  the  high  degree  of  turbulence 
generated  at  separation  overwhelms  all  other  influences  in 
the  approaching  flow. 


d.  MEAN  VELOCITY  AND  TURBULENCE  INTENSITY 

Typically  the  maximum  negative  velocities  in  the  pri¬ 
mary  recirculation  zone  are  10  to  15%  of  the  inlet  center¬ 
line  value  [11,12,15].  This  is  substantially  lower  than  the 
maximum  reverse  flow  velocities  found  in  two-dimensional 
steps.  Eaton  and  Johnston  [6]  quote  values  of  .250^  and 
Bremmer,  et  al .  [3]  quote  a  maximum  reverse  velocity  of 


.  20U^  The  difference  in  these  values  may  be  due  to  step 
height  differences  in  the  experiments. 

Centerline  velocity  decay  following  an  abrupt  axisym- 
metr ic  expansion  is  similar  to  that  for  a  free  jet  in  the 
"core"  region  of  the  inlet  flow.  This  "cote"  region  is 
approximately  five  step  heights  long  for  a  sudden  expansion 
with  an  area  ratio  of  2.  Beyond  this  point  the  similarity 
between  the  free  jet  and  sudden  expansion  ceases  to  exist, 
since  the  free  jet  velocity  decreases  to  zero  while  the  sud¬ 
den  expansion  flow  decays  to  a  fully  developed  turbulent 
pipe  flow.  The  centerline  velocity  decay  appears  to  be 
nearly  independent  of  area  ratio  and  inlet  velocity  profile 
when  plotted  against  non-dimensional  step  height  [15J. 

The  centerline  local  turbulence  intensity  rapidly 
increases  from  a  low  value  (typically  1-2%)  at  the  sudden 
expansion  to  a  value  near  40%  approximately  13  step  heights 
downstream.  Maximum  turbulence  intensities  in  the  shear 
layer,  normalized  with  respect  to  the  inlet  centerline  velo¬ 
city,  are  typically  about.  25%  [3,6,15).  The  reattachment 
length  in  an  ax isymmetr ic  sudden  expansion  is  normally 
between  7  and  9  step  heights  and  appears  to  be  independent 
of  area  ratio  and  inlet  conditions  for  turbulent  flows. 


ANALYTICAL  MODELING  OF  AX  I SYMMETRIC  SUDDEN  EXPANSION 
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Turbulent  sudden  expansion  flows  have  received  varied 
analytical  treatment  in  recent  years.  Probably  the  most 

popular  computer  codes  for  this  turbulent  flow  problem  util-  ; 

lze  the  two-equation  k-e  model  originally  developed  by  Har¬ 
low  and  Nakayma  [19]  which  has  been  modified  by  Launder  and  » 

Spalding  [20]  and  Launder,  et.  al  [21].  Both  Moon  and  Rud- 
mger  [12]  and  Stevenson,  et  al.  [lb]  generated  predictions 
for  axisymetric  sudden  expansions  using  codes  of  this  type. 

The  two  dimensional,  time  averaged  conservation  equations  in 
elliptic  form  were  solved  using  refined  finite  difference 
techniques.  The  two-equation  turbulence  model  requires 
"universal"  turbulence  coefficients.  Moon  and  Rudinger 
demonstrated  that  these  coefficients  were  not,  in  fact, 
universal,  at  least  not  for  recirculating  flows.  Stevenson, 
et  al.  [15]  found  that  although  the  k-€  model  was  adequate 
for  engineering  purposes,  it  did  not  yield  a  precise 
representation  of  the  flow  field.  Gosman,  Khalil  and  Whi- 
telaw  [22]  felt  that  the  dissipation  equation  caused  at 
least  part  of  the  deficiency  in  the  model.  More  complex 
models  (Reynolds  stress),  however,  apparently  do  not  result 
in  any  better  representation  of  the  flow  field  [22]. 
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3.  REACTING  PLOWS 


a.  INTRODUCTION 

The  interaction  between  combustion  and  turbulence  has 
been  the  subject  of  much  research  in  the  past  few  years. 
The  effects  of  combustion  on  the  structure  of  turbulence  and 
the  effects  of  turbulence  on  combustion  (i.e.  chemical  reac¬ 
tion  rates)  is  of  great  interest  to  the  combustion  engineer. 
Knowledge  of  these  effects  permit  the  design  of  more  energy 
efficient  combustion  systems  with  low  pollutant  emission 
while  minimizing  costly  test  programs. 

The  objective  of  most  recent  experimental  work  in  tur¬ 
bulent  reacting  flows  is  to  obtain  data  for  comparison  with 
the  numerical  prediction  codes  that  model  turbulence, 
combustion  and  heat  transfer.  LDV  measurements  in  reacting 
flows  have  been  made  in  various  geometries.  The  more  common 
include:  1)  co-axial  jets,  2)  diffusion  flame  jets,  3) 
bluff  body  flame  holders  in  ducts,  4)  industrial  furnaces, 
and  5)  two-dimensional  rearward  facing  steps.  All  of  these 
geometries  have  a  common  feature  in  that  they  induce  a 
separated  flow  region  which  acts  as  an  anchor  for  the  flame 
front.  High  turbulence  intensities  and  recirculation  zones 
are  characteristic  features  of  these  flows.  Comparisons  of 
mean  velocities,  turbulence  intensities  and  recirculation 
zone  sizes  have  been  made  in  the  above  mentioned  flow  fields 
with  and  without  combustion.  These  comparisons  have  been 


made  in  an  effort  to  understand  the  turbulence-combustion 
interaction.  Turbulence  structure  parameters  (i.e.  skewness 
and  flatness)  are  also  often  compared  to  give  some  insight 
to  the  interaction  process.  Results  from  these  studies  will 
be  discussed  in  the  following  section.  The  effects  of  axial 
pressure  gradient  on  the  turbulent  structure  of  reacting 
flows  will  also  be  discussed.  Pinally,  a  review  of  the 
current  computer  codes  used  to  predict  turbulent  flow  with 
combustion  will  be  presented. 

b.  COMPARISON  OF  FLOWS  WITH  AND  WITHOUT  COMBUSTION 

A  comparative  study  of  cold  and  reacting  flow  around  a 
bluff  body  flame  stabilizer  in  a  duct  was  made  by  Fujii  and 
Eguchi  [24].  Mean  streamwise  velocities  and  turbulence 
intensities  were  measured  downstream  of  the  stabilizer.  A 
homogeneous  propane-air  mixture  was  used  in  the  reacting 
flow  study.  They  found  the  recirculation  zone  to  be  50% 
longer  and  30%  wider  in  the  flow  with  combustion  as  compared 
to  the  isothermal  flow.  The  recirculation  zone  varied  from 
5.4  to  7.5  bluff  body  half-widths  downstream  of  the  stabil¬ 
izer.  This  recirculation  zone  length  was  found  to  be  a 
function  of  the  equivalence  ratio,  attaining  a  minimum  value 
at  stoichimetr ic  fuel-air  ratios  and  reaching  a  maximum 
value  at  the  weak  and  rich  limits.  Local  turbulence  levels 
in  the  reacting  flow  were  found  to  be  much  lower  in  the 
reacting  flow,  typically  50%  of  the  values  found  in  the 
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isothermal  flow.  Fuju  and  Eguchi  proposed  that  this 
suppression  of  turbulence  by  combustion  resulted  from  dila¬ 
tation  by  heat  release  competing  with  turbulence  energy  pro¬ 
duction  in  the  shear  layer.  The  structure  parameters,  skew¬ 
ness  and  flatness,  were  found  to  be  highly  distorted  in  the 
reacting  flow  indicating  highly  anisotropic  turbulence. 
Using  spectrum  analysis  they  also  found  that  the  distinct 
eddy  formation  and  shedding  mechanism,  characteristic  of 
isothermal  separated  flows,  disappeared  in  the  reacting 
flow.  This  was  also  observed  previously  by  Williams,  et  al. 
(25).  The  disappearance  of  eddy  shedding  in  flows  with 
combustion  can  be  explained  by  considering  the  static  pres¬ 
sure  distribution  downstream  of  the  bluff  body  (flame 
holder).  The  minimum  static  pressure  occurs  at  a  point 
immediately  downstream  of  separation.  The  static  pressure 
then  starts  to  rise  as  you  proceed  downstream  in  the  wake 
region  until  a  maximum  is  reached  (where  minimum  gas  velo¬ 
city  occurs).  The  distance  between  separation  and  the  loca¬ 
tion  of  maximum  static  pressure  approximately  defines  the 
eddy  region  length.  With  combustion,  the  eddy  shedding  pro¬ 
cess  is  stopped  because  the  position  of  maximum  static  pres¬ 
sure  shifts  (due  to  the  f undamental  pressure  loss  associated 
with  heat  release)  considerably  closer  to  the  bluff  body  as 
compared  to  the  isothermal  flow. 

Durao  and  Whitelaw  [26]  measured  the  centerline  axial 
mean  velocity  and  corresponding  normal  stress  in  both  disc 


stabilized  diffusion  and  premixed  flames  using  an  LDV .  They 
found  that  the  recirculation  zone  lengths  in  the  isothermal 
flow,  diffusion  flame  and  premixed  flame  cases  were  similar. 
This  is  in  disagreement  with  the  results  of  Pujii  and  Eguchi 
and  may  be  due  to  the  fact  that  these  flames  were  uncon- 
fined.  This  demonstrates  the  need  for  research  to  determine 
the  effect  of  confinement  on  flames.  The  maximum  reverse 
velocity  in  the  recirculation  zone  was  found  to  depend  on 
initial  velocity  and  reached  a  value  approximately  40%  of 
the  inlet  velocity  for  all  flow  cases.  The  normal  stresses 
observed  in  the  attached  flames  were  generally  found  to  be 
lower  than  in  the  isothermal  jet.  This  suggests  that  the 
heat  from  the  surrounding  flame  suppresses  turbulence 
energy.  The  velocity  probability  distribution  functions 
were  found  to  be  near  Gaussian  in  shape  along  the  flame 
centerline.  This  result  would  be  expected  due  to  the  sym¬ 
metry  of  the  flow  and  the  fact  that  the  large  gradients  and 
flow  intermittenc les  occur  at  the  f lame  boundary. 

In  a  study  of  enclosed  turbulent  diffusion  flames, 
Hartmann  [27]  made  LDV  measurements  in  two  different  combus¬ 
tion  chambers.  In  both  chambers  either  propane  or  natural 
gas  was  introduced  through  a  central  orifice  surrounded  by 
ait  passing  through  an  annular  swirler.  The  mean  axial 
velocities  and  rms  fluctuations  were  measured  using  an  LDV 
employing  frequency  shifting  operating  in  the  backscatter 
mode.  The  backscatter  mode  was  used  for  two  reasons:  1) 
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optical  access  was  available  on  one  side  of  the  chamber  only 
and  2)  the  backscatter  mode  was  believed  to  minimize  the 
problems  of  beam  bending  and  optical  alignment  due  to 
refractive  index  gradients  in  the  flow  field.  Hartmann 
found  that  mean  streamwise  velocities  increased  in  the  flows 
with  combustion,  as  expected,  due  to  the  reduction  in  gas 
density.  The  recirculation  zone  was  found  to  be  narrower 
(approximately  70 %  of  the  value  found  in  the  isothermal 
case)  in  the  flow  with  combustion.  The  local  turbulence 
intensity  was  typically  found  to  be  15%  lower  in  the  com¬ 
busting  flow  as  compared  to  the  isothermal  flow.  This  is  in 
agreement  with  both  Fu^ii  and  Eguchi  [24]  and  Durao  and  Whi- 
telaw  [26]. 

Baker ,  et  al.  [28]  measured  the  three  components  of 
mean  velocity  and  corresponding  normal  stresses  in  an 
axi3ymmetric  coaxial  jet  fu;  ,ace  with  and  without  combustion 
using  a  one  component  LDV.  Comparisons  of  the  flow  field 
with  and  without  swirl  were  also  made.  Profiles  of  the 
measured  turbulent  kinetic  energy  in  the  isothermal  flows 
were  presented  and  the  regions  of  near - isotropic  turbulence 
were  identified.  The  study  showed  that  the  recirculation 
regions  were  substantially  different  in  size  for  combusting 
flow  as  compared  to  isothermal  flow  and  that  the  turbulence 
was  significantly  more  anisotropic  over  most  of  the  flow 
field  as  evidenced  by  the  highly  skewed  velocity  probability 
distribution  functions.  The  recirculation  zone  in  this 
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study  was  found  to  increase  in  length.  This  is  contrary  to 

the  results  of  Fujii  and  Eguchi  [24]  when  considering  that 

/ 

the  recirculation  zones  are  toward  the  outside  of  the 
combustor  in  the  present  study  (similar  to  an  axisymmetric 
sudden  expansion)  but  toward  the  center  of  the  test  section 
in  the  flow  around  a  bluff  body  flame  stabilizer.  Integrat¬ 
ing  the  rms  velocity  fluctuations  across  the  radius  of  the 
furnace  in  both  the  isothermal  and  reacting  flows  indicated 
that  the  velocity  fluctuations  (turbulent  kinetic  energy) 
increased  significantly  (125%)  as  a  consequence  of  the 
combustion.  The  explanation  for  this  result,  which  is  in 
disagreement  with  the  four  previously  discussed  studies,  is 
not  clear.  The  technique  of  integrating  the  rms  fluctua¬ 
tions  across  the  radius  and  directly  comparing  the  two  flow 
fields  may  be  questionable,  however.  It  would  seen  that  the 
energy  added  from  the  heat  input  in  the  reacting  case  would 
have  to  be  included  in  the  turbulent  kinetic  energy  balance 
in  order  to  compare  the  two  flow  fields  directly. 

In  an  effort  to  evaluate  the  effect  of  combustion  on 
the  turbulent  structure  of  a  two-dimensional  rearward  facing 
step  flow,  Pitz  {29]  made  extensive  LDV  measurements  in  such 
a  flow  with  and  without  combustion.  Detailed  mappings  of 
mean  streamwise  velocities  and  turbulence  intensities  were 
presented.  The  large  scale  turbulence  structure  was  observed 
using  high  speed  Schlieren  photography  and  spectral  analysis 
(LDV).  A  homogeneous  propane-air  mixture  with  an 
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equivalence  ratio  of  0.57  was  used  in  the  combustion  study. 
The  reattachment  length  in  the  flow  with  combustion  was 
found  to  be  30%  shorter  than  that  found  in  the  isothermal 
flow.  The  recirculation  zone  was  also  found  to  be  thinner; 
that  is,  the  shear  layer  boundary  shifted  towards  the  wall 
in  the  reacting  flow.  This  was  due  to  the  increased  velo¬ 
city  in  the  core  flow  at  the  top  of  the  shear  layer  (due  to 
combustion)  which  caused  the  upper  boundary,  as  defined  by 
the  location  of  zero  mean  velocity,  to  shift  toward  the 
wall.  Peaks  in  the  turbulence  intensity  profiles  confirmed 
this  shift  of  the  recirculation  zone  boundaries.  It  is  of 
interest  to  note,  however,  that  Schlieren  visualization 
showed  the  f lame  boundary  propagating  further  into  the  core 
region  than  the  shear  layer  boundary  (defined  by  the  mean 
velocity  profiles).  The  maximum  normalized  turbulence  level 
in  the  reacting  flow  was  found  to  be  approximately  30% 
higher  than  that  in  the  isothermal  flow  and  was  located  in 
the  shear  layer  one  step  height  downstream  of  the  step. 
Although  the  maximum  normalized  turbulence  level  was  higher 
m  the  reacting  flow  case,  local  turbulence  intensities  were 
found  to  be  lower  due  to  the  higher  local  velocities  present 
in  the  reacting  flow.  Typically,  the  values  of  local  tur¬ 
bulence  intensity  throughout  the  flow  were  a  factor  of  three 
less  in  the  reacting  case  as  compared  to  the  isothermal  case 
at  a  plane  in  the  vicinity  of  reattachment.  This  seems  to 
suggest  that  more  turbulence  is  generated  in  the  shear  layer 
in  the  reacting  flow  case,  but  is  suppressed  downstream  due 
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to  the  combustion  process.  The  reason  for  increased  shear 
generation  in  the  reacting  flow  case  is  not  clear.  In 
agreement  with  Fujii  and  Eguchi  [24]  and  Williams,  et  al. 
[25],  combustion  was  found  to  reduce  the  eddy  coalescence  in 
the  reacting  shear  layer.  The  pairing  process  was  nearly 
eliminated  as  a  growth  mechanism  in  the  reacting  layer. 

c.  EFFECTS  OF  AXIAL  PRESSURE  GRADIENTS 

It  was  mentioned  earlier  that  axial  pressure  gradients 
ate  suspected  of  playing  a  large  role  m  determining  the 
location  of  reattachment  in  isothermal  ax  1 symmetr ic  sudden 
expansions  (4J.  The  effects  of  axial  pressure  gradients  on 
turbulent  diffusion  flames  were  studied  by  Starner  and 
Bilger  [30],  In  this  study,  various  positive  and  negative 
pressure  gradients  were  imposed  on  a  coaxial  jet  configura¬ 
tion  and  LDV  measurements  were  made.  The  measurements 
included  axial  and  radial  mean  velocities  and  turbulence 
intensities.  Reynolds  stress  correlations  along  with  skew¬ 
ness  and  flatness  parameters  were  also  presented.  It  should 
be  noted  that  the  turbulence  structure  in  confined  coaxial 
jet  flows  is  much  more  complex  than  that  of  the  flow  in  an 
axisymmetric  sudden  expansion.  However,  the  primary  to 
secondary  jet  velocity  ratio  used  in  this  experiment 
(Up/ug  »  10)  causes  the  flow  structure  to  approach  that  of  a 
sudden  expansion  as  opposed  to  an  ejector.  They  found  that 
local  centerline  turbulence  intensities  increased  by 
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approximately  a  factor  of  two  in  flows  with  either  positive 
or  negative  pressure  gradients  imposed  on  them  (for  dif¬ 
ferent  reasons)  as  compared  with  the  zero  pressure  gradient 
case.  This  result  is  very  interesting  as  Glass  and  Bilger 
[31]  found  that  the  normalized  turbulence  intensity  was  much 
the  same  for  both  isothermal  and  reacting  flows  with  a  near 
zero  pressure  gradient.  This  indicates  that  the  normal 
shear -gener at ed  turbulence  mechanism  is  dominant  in  both 
isothermal  and  reacting  flows  with  near  zero  axial  pressure 
gradient.  The  increase  in  turbulence  intensity  for  adverse 

pressure  gradient  flows  indicate  that  the  advection  term, 
0  k 

i.e.  pu(-j^),  is  of  the  same  order  as  the  production  term 
(at  the  radius  of  maximum  shear  stress)  in  the  Favre  aver¬ 
aged  turbulent  kinetic  energy  equation  [32].  In  accelerat¬ 
ing  flows  it.  appeared  that  the  turbulence  was  produced  by 
the  combustion  process.  This  was  suggested  as  the  ratio  of 

the  source  (or  sink)  term  of  turbulent  kinetic  energy,  i.e. 
3p 

u  to  the  production  term  (at  radius  of  maximum  shear 
stressi  increased  m  the  downstream  direction.  Purther  work 
in  this  area  is  indicated  as  there  are  significant  implica¬ 
tions  tor  turbulence -combust  ion  modelling.  A  detailed  map¬ 
ping  of  the  complete  flow  field  with  a  two-component  LDV 
system  may  resolve  more  of  the  terms  in  the  turbulent 
kinetic  energy  balance 
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These  observations  could  explain  why  some  experimenters 
found  turbulence  to  increase  with  combustion  while  others 
found  turbulence  to  be  suppressed  by  combustion.  One  should 
be  careful  when  stating  that  a  complete  reacting  flow  falls 
into  one  of  these  two  broad  categories  (  i.e.  suppression  or 
generation  of  turbulence  due  to  the  combustion  process). 
Different  turbulence  mechanisms  (i.e.  production,  dissipa¬ 
tion,  advection,  etc.)  dominate  in  different  regions  of  the 
flow.  For  example,  far  downstream  of  the  step  in  an  ax  1  sym¬ 
metric  sudden  expansion  turbulence  production  is  insignifi¬ 
cant,  while  turbulence  production  is  very  important  in  the 
turbulent  kinetic  energy  balance  immediately  downstream  of 
the  step.  When  comparing  turbulence  character ist ics  of  dif¬ 
ferent  flow  geometries,  one  must  be  sure  the  fluid  mechan¬ 
isms  in  the  two  regions  of  concern  are  sufficiently  similar. 

The  decay  of  centerline  velocity  was  also  found  to 
depend  strongly  on  pressure  gradient.  Typically  the  more 
positive  the  pressure  gradient  (diffusing  flow)  the  faster 
the  decay.  Reacting  flows  were  found  to  be  more  sensitive  to 
imposed  pressure  gradients  due  to  the  low  density  of  the  hot 
gases.  The  structure  parameters  (skewness  and  flatness)  and 
the  correlation  coefficients  were  found  to  be  independent  of 
pressure  gradient.  This  would  imply  that  the  more  compli¬ 
cated  Reynolds  stress  turbulence  models  (that  take  into 
account  the  anisotropic  nature  of  turbulence)  may  not  be 
necessary  ( 30 ] . 


2. 5 


d.  NUMERICAL  MODELING 


One  important  motivation  for  making  detailed  LDV  meas¬ 
urements  in  turbulent  reacting  flows  is  to  generate  experi¬ 
mental  data  which  can  be  used  to  verify  the  predictions 
obtained  from  computer  codes.  Most  measurements  of  this 
type  have  been  aimed  at  obtaining  mean  velocity  and  tur¬ 
bulence  intensity  data  throughout  the  flow  field  as  men¬ 
tioned  in  Section  3.1. 

The  most  popular  numerical  prediction  code  in  use  today 
employs  the  two-equation  turbulence  model  (k  -  e)  to  com¬ 
plete  closure  for  the  time  averaged  Nav ier -Stokes  equations. 
Typically  one  of  three  combustion  models  is  used  character¬ 
ized  by:  1)  instant  reaction,  2)  instant  reaction  with 
scalar  fluctuations,  or  3)  Arrhenius  reaction  or  eddy-break 
up  with  scalar  fluctuations. 

Khalil,  Spalding  and  Whitelaw  (33)  performed  numerical 
predictions  with  a  code  as  described  above  that  used  all 
three  combustion  models  and  obtained  results  that  agreed 
with  available  experimental  mean  velocity  data  [28]  to  a 
good  degree  over  most  of  the  flow  field.  The  maximum 
disagreement  occurred  at  the  centerline  of  the  flow  which  is 
also  where  prediction  problems  occurred  in  isothermal  flows 
fl5).  The  instant  reaction  with  scalar  fluctuations  combus¬ 
tion  model  was  found  to  give  marginally  better  predictions 
than  the  other  two  models.  The  relatively  simple  (k  -  c) 
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turbulence  model,  as  compared  to  the  more  complicated  tur¬ 
bulence  stress  model,  was  found  to  predict  the  flow  suffi¬ 
ciently  well  in  view  of  the  additional  complexities  involved 
in  employing  the  more  accurate  stress  model  [33].  They 
found,  in  fact,  that  anisotropic  turbulence  models  will  not 
be  useful  until  improvements  are  made  in  the  existing 
combust  ion  models.  Additional  numerical  predictions  in  com¬ 
busting  flows  were  made  by  Hutchinson,  Khalil  and  Whitelaw 
(34)  and  Gosman,  Lockwood  and  Salooja  [3b]  with  similar 
results . 


Bray  [32]  has  questioned  the  common  practice  of  using 
constant  density-empirical  closures  and  modeling  equations 
in  reacting  flow  situations.  Use  of  these  models  is  hard  to 
justify  theoretically  when  it  is  realized  that  many  closures 
originated  from  simple  dimensional  analysis,  while  combus¬ 
tion  introduces  additional  dimensionless  groups  such  as  den¬ 
sity  ratio.  These  arguments  are  valid  from  the  theoretical 
standpoint,  but  until  better  combustion  models  are  developed 
constant  density  closure  models  seem  adequate.  The  most 
important  need  of  the  combustion  modeler  at  the  moment  is 
reliable  measurements  of  velocity,  temperature,  species  con¬ 
centration  and  the  corresponding  correlations  in  well 
defined  turbulent  reacting  flows. 
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4. 


LDV  MEASUREMENT  ACCURACY 


a.  INTRODUCTION 

The  LDV  has  many  desirable  qualities,  as  noted  in  Sec¬ 
tion  I,  which  make  it  a  useful  diagnostic  tool  m  highly 
turbulent  flows.  In  reacting  flows  it  is  especially  advan¬ 
tageous  due  to  the  hostile  environment  present  in  these 
flows  which  prohibits  the  use  of  conventional  velocity 
measuring  instruments. 

The  individual  realization  LDV  is,  in  principle,  an 
absolute  instrument;  that  is,  it  requires  no  calibration. 
However  the  finite  size  of  the  probe  volume,  the  fact  that 
the  measured  velocity  is  actually  the  velocity  of  a  small 
particle  as  it  passes  through  the  probe  volume  (as  opposed 
to  the  true  gas  velocity),  and  the  requirements  of  statisti¬ 
cal  averaging,  especially  in  highly  turbulent  flows,  may 
introduce  errors  into  the  measurements.  The  primary  sources 
of  error  introduced  by  this  instrument  are  discussed  m  more 
detail  below.  A  discussion  of  velocity  bias  and  guidelines 
for  the  correction  of  biased  velocity  data  are  also 
pr  esented . 


b.  LDV  MEASUREMENT  ERRORS 

Several  sources  of  error  have  been  discovered  which  are 
inherent  in  LDV  measurements.  These  errors  can  result  from 
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particle  seeding  effects,  probe  volume  effects,  the  data 
sampling  method  used,  and  signal  processor  operating  charac¬ 
teristics.  Many  of  these  errors  have  been  fully  analyzed 
and  are  easily  eliminated  or  reduced  to  an  insignificant 
level.  A  brief  review  of  the  more  significant  errors  is 
presented  below.  The  effect  of  refractive  index  gradient 
effects  in  reacting  flows  is  also  discussed. 

Two  primary  sources  of  error  can  affect  the  accuracy  of 
LDV  measurements  made  m  highly  turbulent  flows.  These  are 
velocity  bias  and  incomplete  signal  bias.  Velocity  bias 
occurs  when  the  velocity  data  are  obtained  at  unequal  time 
intervals  as  controlled  by  random  particle  passages  through 
the  probe  volume.  Since  more  particles  per  unit  time  pass 
through  the  probe  volume  during  time  intervals  when  the 
velocity  is  high,  the  mean  velocity  calculated  by  a  simple 
averaging  of  the  data  is  higher  than  the  true  mean.  Tur¬ 
bulence  intensity  and  the  higher  order  moments  are  also 
affected.  Although  several  correction  schemes  have  been 
proposed  for  the  reduction  of  velocity  bias,  no  one  method 
has  yet  been  selected  as  a  standard.  Sampling  techniques 
have  been  proposed  that  eliminate  velocity  bias,  but  these 
techniques  require  high  seeding  densities  which  may  not  be 
achievable  in  some  flows.  Since  the  velocity  bias  error  can 
be  significant  in  flows  of  the  type  under  consideration 
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here,  it  is  important  to  account  for  it.  A  brief  review  of 

% 

the  status  of  velocity  bias  will  be  presented  in  the  follow¬ 
ing  section. 

Incomplete  signal  bias  occurs  when  particles  traverse 
the  fringes  of  the  probe  volume  at  large  angles  of  attack 
and  generate  an  insufficient  number  of  fringe  crossings  to 
allow  signal  validation  by  the  processor.  This  also  leads 
to  an  incorrectly  high  mean  velocity.  A  detailed  discussion 
of  other  bias  errors  which  may  occur  in  LDV  measurements  is 
given  by  Thompson  and  Flack  [2j.  These  errors  include  fre¬ 
quency  broadening,  directional  ambiguity,  clock  errors  in 
high  speed  counters,  particle  distribution  bias  resulting 
from  non-uniform  seed  distribution,  particle  lag  bias  and 
particle  acceleration  bias.  There  can  also  be  a  bias  if 
there  is  a  steep  velocity  gradient  in  the  probe  volume. 
Usually  these  errors  are  small.  For  the  measurements  of 
interest  here,  only  incomplete  signal  bias  and  velocity  bias 
need  to  be  considered  and,  as  shown  in  Ref.  [4]  incomplete 
signal  bias  is  easily  eliminated  by  appropriate  frequency 
shifting. 

Refractive  index  gradients  generated  by  the  large  den¬ 
sity,  temperature  and  species  concentration  gradients  in 
reacting  flows  can  be  a  concern  when  using  an  LDV  system. 
These  gradients  not  only  cause  the  incoming  laser  beams  to 
bend  due  to  refraction,  but  also  cause  the  beams  to  diverge 
and  lose  spatial  coherence  giving  rise  to  optical  alignment 
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problems  and  poor  s ignal-to-no lse  ratio.  To  date,  no  cri¬ 
teria  exist  to  analytically  predict  whether  refractive  index 
gradients  will  cause  problems.  However,  there  is  empirical 
evidence  on  the  length  of  a  combustion  chamber  that  the 
beams  can  pass  through  and  still  yield  good  signals. 

In  a  large  industrial  furnace  (2  m.  square)  Baker, 
Hutchinson  and  Whitelaw  (36]  measured  velocities  and  tur¬ 
bulence  intensities  using  an  LDV.  They  found  that  the  beam 
waist  diameter  increased  by  about  a  factor  of  four  as  the 
beams  passed  through  the  furnace.  The  beams  also  fluctuated 
approximately  2.5  mm  about  a  .jean  position  due  to  the  fluc¬ 
tuating  refractive  index  in  the  flow.  Therefore  the  beams 
crossed  intermittently,  causing  severe  signal  distortion. 
This  is,  in  fact,  the  only  published  literature  that  identi¬ 
fies  refractive  index  gradients  as  being  a  serious  problem. 

Typically  the  signals  are  of  good  enough  quality  to 
give  accurate  TJDV  measurement,  especially  in  small  diameter 
furnaces  [29].  Careful  alignment  of  the  probe  volume  image 
on  the  pinhole  of  the  photomultiplier  tube  after  the  beams 
have  passed  through  the  combusting  flow  field  will  ensure 
good  signal  quality.  Operating  in  the  backscatter  mode  has 
been  suggested  (27]  to  minimize  the  path  length  traveled  by 
the  beams.  Another  advantage  of  this  method  is  that  the 
scattered  light  is  transmitted  back  along  approximately  the 


same  path  the  incoming  beams  traverse.  This  means  that 
refraction  occur  mg  as  the  beams  go  into  the  flow  happens  in 
the  reverse  order  as  the  signal  comes  out. 

As  a  final  note,  Barlow  [37]  demonstrated  in  a  1.2 
meter  diameter  furnace  that  photon  correlation  methods  are 
well  suited  to  measure  signals  that  are  intermittent  or 
weak.  Prom  these  and  other  experimental  observations  it 
appears  that  the  refraction  problem  is  not  serious  in 
combustion  studies  near  atmospheric  pressure  for  test  zones 
of  0.5  meters  or  less  and  that  measurments,  albeit  of  lower 
quality,  can  be  made  in  flows  up  to  about  1.5  meters  in 
cross-section. 

C.  VELOCITY  BIAS 

McLaughlin  and  Tiederman  [38]  proposed  in  1973  that 
measurements  with  an  individual  realization  LDV  would  con¬ 
tain  a  bia3  commonly  referred  to  as  velocity  or  sampling 
bias.  They  proposed  a  weighting  factor  correction  using  the 
reciprocal  of  the  velocity  components.  Since  that  time 
there  have  been  several  analytical  studies  with  a  variety  of 
proposed  corrections  for  velocity  bias.  Those  are  summar¬ 
ized  in  Table  2.  Some  recent  experimental  work  related  to 
the  velocity  bias  problem  is  summarized  in  Table  3.  The 
tables  indicate  that  researchers  use  four  characteristic 
time  scales  in  an  effort  to  quantify  the  biasing  effect. 
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TABLE  2.  THEORETICAL  REVIEW  OF  VELOCITY  BIAS 
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than  <t  >  measurement  every  sample  period 


ABLE  3.  EXPERIMENTAL  REVIEW  OF  VELOCITY  BIAS 


2  )  the 


They  are:  1)  the  time  between  particle  arrivals,  r  , 
residence  time  of  the  particles  in  the  probe  volume,  r  ,  3) 

the  sampling  time,  t  ,  and  4)  the  characteristic  time  scale 
of  turbulence,  r  . 

Attempts  to  experimentally  verify  that  velocity  bias 
does  in  fact  exist  were  largely  unsuccessful  until  the  work 
of  Roesler,  et  al.  [4]  in  1980  in  which  the  existence  of 
both  incomplete  signal  bias  and  velocity  bias  were  clearly 
shown  experimentally.  More  recently  Johnson,  et  al.  [46] 
have  also  experimentally  verified  the  existence  of  velocity 
bias,  and  have  proposed  a  two  dimensional  velocity  weighting 
factor  to  correct  it.  They  also  tried  using  a  residence 
time  weighting  factor  but  had  little  success.  No  reason  was 
given  for  the  failure  of  this  technique,  but  the  data  pro¬ 
cessing  is  suspected  as  the  problem.  The  two-dimensional 
velocity  weighting  factor  appears  to  be  the  best  method  of 
correcting  for  velocity  bias  when  2- D  measurements  are  avail¬ 
able  and  unbiased  data  cannot  be  obtained  using  the  equal 
time  interval  sampling  method  of  Roesler,  et  al.  [4J. 

In  an  effort  to  determine  the  potential  errors  induced 
by  biased  sampling,  Tiederman  [47]  applied  both  one- 
dimensional  and  two-dimensional  correction  factors  to  biased 
two-d imens ional  LDV  data  obtained  from  the  IT  wind  tunnel  at 
Arnold  Engineering  Development.  Center.  Figures  2  and  3  show 
the  biased  (c*>  =  1)  and  corrected  mean  velocity  and  rms  fluc¬ 
tuating  velocity  profiles  respectively  in  a  separated  region 


F inure  2.  Corrected  and  Uncorrected  (lean 

Streamwise  Velocity  Profiles  in  a 
Separated  Region  (from  Ref.  [47]) 


Figure  3.  Corrected  and  Uncorrected  rms  Streanwise 
Velocity  Profiles  in  a  Separated  Region 
( from  Ref .  [47]) . 
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downstream  of  a  hump.  The  trends  of  the  mean  velocity  data 
(biased  vs.  2-D  correction)  are  as  expected.  The  1-D 
correction  scheme  was  found  to  give  spurious  results  due  to 
the  singularities  that  occurred  when  the  velocity  component 
equaled  zero.  These  singularities  were  found  to  be  absent 
when  using  the  2-D  correction  scheme.  The  profiles  (biased 
vs.  2-D  correction)  of  rms  fluctuating  velocity  were  found 
to  be  s ign if  leant ly  different  as  shown  in  Pigure  3.  This 
result  is  important  as  the  spatial  distribution  of  tur¬ 
bulence  intensity  is  altered  by  a  sizable  amount  when  biased 
velocity  data  is  used.  It  appears,  however,  that  the  max¬ 
imum  rms  fluctuation  velocity  is  not  significantly  changed. 


d.  BIAS  CORRECTION 


Based  on  the  experimental  and  theoretical  findings  to 
date,  a  set  of  guidelines  for  velocity  and  incomplete  signal 
bias  elimination  or  correction  can  be  formulated. 


a.  Velocity  and  incomplete  signal  biases  are  only 
important  in  highly  turbulent  flows.  Therefore, 
no  correction  is  needed  for  flows  with  local 
turbulence  intensity  levels  below  about  10%. 

b.  Incomplete  signal  bias  can  be  virtually  eliminated 
by  frequency  shifting,  and  since  frequency  shifting 
is  essential  in  highly  turbulent  flows  to  resolve 
the  flow  direction,  that  bias  is  of  little  concern 
in  most  appl icat ions . 

c.  The  easiest  way  to  eliminate  velocity  bias  is  to 
heavily  seed  the  flow  and  control  the  processor 
sampling  rate  such  that  nearly  time  averaged  data 
are  collected.  This  sampling  technique  is  recom¬ 
mended  whenever  it  is  possible  to  use  it.  The 
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seeding  density  should  provide  a  particle  arrival 
rate  two  orders  of  magnitude  or  more  greater  than 
the  sampling  rate.  The  technique  is  described  in 
Ref.  (4]. 

d.  Velocity  bias  corrections  of  data  at  local  turbu¬ 
lence  intensity  levels  up  to  20%  should  be  made 
only  if  high  levels  of  accuracy  are  required.  For 
flows  in  which  the  mean  velocity  is  not  near  zero 
and  the  local  turbulence  is  of  the  order  of  25%  or 
Less,  the  one-dimensional  McLaugh 1 m-T leder man 
correction  appears  to  work  well  for  the  main  veloc¬ 
ity  component.  [f  two-component  measurements  are 
made  1  lie  two-dimensional  correction  scheme  should 
be  user], 

e.  Unless  velocity  bias  is  eliminated  by  equal  time 
interval  sampling  (at  high  seed  density)  or  the 
velocity  sampling  is  controlled  by  particle 
arrival  statistics  such  that  "completely  biased" 
data  is  obtained,  the  amount  of  bias  in  a  data 
set  is  unknown.  Obviously  it  is  inappropriate  to 
use  any  of  the  proposed  bias  correction  methods 
on  a  "partially  biased"  data  set,  since  the 
resulting  error  may  exceed  that  in  the  original 
data.  Unfortunately  much  of  the  data  in  the 
published  literature  is  ill-defined  in  this  sense 
and  must  be  used  with  this  limitation  in  mind. 
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SECTION  III 


1 


EXPERIMENTAL  APPARATUS 


1.  INTRODUCTION 

The  experimental  mapping  of  the  flow  field  in  the 
axisyrametric  sudden  expansion  required  various  types  of 
instrumentation  and  hardware.  The  LDV  used  in  this  experi¬ 
ment  was  designed  specifically  to  allow  investigation  of  the 
effect  of  data  acquisition  methods  and  optical  parameters  on 
measurements  in  highly  turbulent  flows.  The  system  allowed 
for  variation  in  LDV  optical  parameters,  seeding  particle 
arrival  rate  and  data  sampling  conditions.  This  section 
describes  the  expr imental  apparatus  which  may  be  divided 
into  five  major  subsystems : 

1.  the  LDV  optical  system 

2.  the  flow  system 

).  the  fuel  system 

4.  the  data  collection,  storage  and  processing 

system 

5.  the  seeding  systems 


40 


J 


THE  LDV  OPTICAL  SYSTEM 


> 


The  LDV  used  is  a  one  component  system  operating  in  the 
dual-beam  forward  scatter  mode.  It  has  the  capability  of 
changing  probe  volume  size,  fringe  spacing,  and  angular 
orientation.  Frequency  shifting  one  or  both  beams  and 
traversing  the  probe  volume  in  three-dimensional  space  are 
also  possible.  A  schematic  drawing  of  the  general  layout  of 
the  LDV  system  is  shown  in  Figure  4. 

Laser  light  fot  the  system  is  provided  by  a  five  watt 
Coherent  Radiation  Model  52  argon  ion  laser  normally 
operated  on  the  green  line  (0.5145  Mm).  (The  laser  power 
was  set  between  100  and  200  raW  for  all  experimental  runs  in 
the  present  study. )  The  beam  exits  the  laser  and  enters  a 
polarization  rotator  (Spectra-Phys ics ,  Model  310-21). 
Rotating  the  plane  of  polarization  of  the  beam  perpendicular 
to  the  beam  dividing  prism  insures  that  the  beam  is  split 
into  two  equal  intensity  beams.  This  produces  maximum 
fringe  contrast  in  the  probe  volume.  Upon  leaving  the 
polarization  rotator,  the  beam  passes  through  a  beam 
expander  telescope  composed  of  a  44  mm  lens  and  a  66  mm 
lens.  Traversing  the  second  lens  over  a  7.5  mm  range  varies 
the  beam  waist  diameter  from  60  to  500  fun. 

Following  the  telescope,  two  broadband  all-dielectr ic 
mirrors  (Newport  Research  Corp. )  direct  the  beam  to  the  beam 
splitter  (TSI,  Model  1916-1)  on  the  upper  table  of  the 
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UPPER  OPTICS  PACKAGE 
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optics  package.  The  beam  splitter  divides  the  entering  beam 
into  two  parallel  equal  intensity  beams  which  enter  two 
acousto-optic  modulators  (Intra-Action  Corp.,  Model  #ADM- 
40).  The  modulators  shift  the  frequency  of  the  incoming 
beam  by  an  amount  equal  to  the  frequency  of  the  driver.  The 
frequency  shift  may  be  either  up  or  down  allowing  for  a  wide 
range  of  net  frequency  shifts  between  the  two  beams.  A  net 
frequency  shift  of  10  MHz  was  used  in  the  present  study. 
This  allowed  unambiguous  measurement  of  negative  velocities 
up  to  approximately  40  ra/s. 

Upon  leaving  the  modulators,  the  beams  are  reflected  by 
adjustable  mirrors  (Newport  Research  Corp.,  Model  #600-2)  to 
a  sliding  prism.  Adjustment  of  the  prism  changes  the  beam 
separation  and  therefore  the  converging  beam  angle,  thus 
controlling  the  fringe  spacing  and  the  number  of  fringes  in 
the  probe  volume.  The  adjustable  mirrors  are  used  to  posi¬ 
tion  the  beams  such  that  they  cross  at  their  waists  on  the 
optical  axis  after  passing  through  the  transmitting  lens. 
The  transmitting  lens  (TSI,  Model  #918)  has  a  focal  length 
of  250  mm. 

Scattered  light  from  particles  passing  through  the 
probe  volume  is  collected,  collimated  and  focused  by  a  pair 
of  receiving  lenses  (TSI,  Models  #917  and  #918)  mounted 
several  centimeters  apart.  The  focal  lengths  of  the  receiv¬ 
ing  lenses  are  250  nan  and  120.6  mm,  r espect ively .  The 
focused  light  is  reflected  by  a  mirror  mounted  on  an 
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adjustable  fixture  (Newport  Reseach  Corp. ,  Model  1600-2). 
This  allows  fine  lateral  adjustment  of  the  focused  spot, 
insuring  that  it  is  located  on  the  200  <im  pinhole.  The 
receiving  optics  package  may  be  moved  along  the  optical  axis 
for  coarse  adjustment  of  the  probe  volume  image  on  the 
pinhole  with  fine  adjustment  provided  by  the  threaded 
pinhole  mounting.  Beam  stops  on  the  initial  receiving  lens 
block  the  direct  laser  beams  and  allow  only  scattered  light 
to  pass  into  the  PM  tube. 

The  upper  transmitting  optics  table  is  mounted  on  bear¬ 
ings  which  allow  angular  rotation  about  the  optical  axis  and 
thus  permit  any  velocity  component  in  a  perpendicular  plane 
to  be  measured.  The  entire  optics  package  (including  laser) 
is  mounted  on  a  3-axis  milling  machine  table.  Three  Bodine 
DC  gearmotors  with  variable  speed  control  are  used  to  drive 
the  mill  table.  Linear  potentiometers  (New  England  Instru¬ 
ments)  with  a  linearity  of  0.25%  are  used  to  obtain  an 
electrical  readout  of  position  on  digital  panel  meters  which 
read  directly  in  millimeters  to  an  accuracy  of  ±  0.1  mm. 
The  traverse  range  is  254  mm  (10  in.)  in  the  vertical  (y) 
direction  and  152  mm  (6  in.)  in  the  x  and  z  directions.  A 
more  detailed  description  of  the  entire  optical  system 
including  the  individual  components  is  given  by  McVey  [48]. 
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THE  PLOW  SYSTEM 


3  . 


The  flow  system  was  designed  to  provide  a  flexible  sys¬ 
tem  allowing  easy  optical  access,  while  fitting  within  the 
diagnostic  range  of  the  LDV  for  both  isothermal  and  reacting 
flows.  The  geometry  used  in  this  study  is  that  of  an 
axisymmetr ic  sudden  expansion.  One  unigue  feature  of  the 
test  rig  is  that  a  cylindrical  quartz  tube  was  used  as  the 
test  section.  This  eliminated  the  flow  field  perturbation 
induced  by  inserting  flat  windows  in  a  cylindrical  pipe. 
Mote  will  be  said  about  the  te3t  section  later.  The  flow 
system  consisted  of  eight  major  part3,  as  shown  in  Pigure  S. 

1.  a  radial  vane  blower 

2.  a  flow  conditioning  section 

3.  a  connecting  duct 

4.  main  fuel  injection 

5.  torch  ignitor 

6.  a  cylindrical  quartz  test  section 

7.  a  thermocouple  rake 

8.  an  extension  duct 

The  radial  vane  blower  was  a  Peerless  Model  PWB4GA 
driven  by  a  variable  speed  direct  current  motor.  The 
blower -motor  combination  allows  a  flow  capacity  of  1100  cfm. 

The  flow  conditioning  section  consisted  of  a  set  of 
flow  str aighteners  76.2  non  (3  in.)  in  diameter.  The  ele¬ 
ments  of  this  section  included  wire  window  screen  followed 
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RADIAL 

BLOWER 


by  a  honeycomb  of  6.35  ram  (0.25  in.)  diameter  plastic  soda 
straws  25.4  mm  (1  in.)  in  length.  This  section  exited  into 
a  series  of  four  window  screens  spaced  25.4  ram  (1  in.) 
apart.  The  conditioning  section  was  178  ram  ( 7  in.)  m 
length  and  was  connected  to  the  blower  via  a  convergent 
adapter . 

The  connecting  duct,  fabricated  from  76.2  mm  (3  in.) 
diameter  standard  schedule  40  steel  pipe,  was  3.032  m  (.  10 
ft.)  in  length.  This  length  was  chosen  to  give  a  fully 
developed  velocity  profile  as  an  inlet  boundary  condition  to 
the  ax lsymmetr ic  sudden  expansion. 

The  main  fuel  (gaseous  propane)  was  in3ected  through 
four  plain  hole  orifices,  located  90°  apart  circumferen¬ 
tially,  in  a  plane  located  2.804  m  (110  in.)  upstream  of  the 
sudden  expansion.  This  allowed  sufficient  mixing  for  flame 
stabilization  in  the  dump-section.  A  completely  premixed 
fuel-air  mixture  was  not  desired  at  the  dump  plane  as  this 
mixture  would  lie  outside  the  burning  limits  for  the  lean 
overall  fuel-air  mixture  used  in  this  experiment.  Using  a 
lean  fuel-air  mixture  allowed  for  steady  state  operation  by 
keeping  the  test  section  walls  below  the  annealing  tempera¬ 
ture  for  quartz  (1070°C). 

A  pitot  tube  was  inserted  in  the  connecting  duct  695  mm 
(:  27  in.)  upstream  of  the  sudden  expansion.  This  was  used 
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to  set  and  monitor  the  inlet  centerline  flow  velocity.  The 
pitot  tube  remained  inserted  in  the  duct  during  all  test 
runs . 


A  propane-air  torch  ignitor  was  used  to  ignite  the 
fuel-air  mixture.  The  ignitor  was  mounted  flush  to  the  con¬ 
necting  duct  wall  to  minimize  flow  field  disturbances  and 
was  located  175  mm  ( :  7  in. )  upstream  of  the  sudden  expan¬ 
sion.  The  torch  ignitor  was  fabricated  from  stainless  steel 
tubing  and  is  shown  in  Pigute  6.  The  propane  flows  down  the 
annulus  formed  by  the  3/8"  tube  and  the  1/4"  ceramic  insula¬ 
tor.  The  air  flows  down  the  annulus  formed  by  the  1/2*  tube 
and  the  3/8*  tube.  The  3/8"  tube  is  approximately  1  inch 
shorter  than  the  1/2"  outer  tube  in  order  to  allow  the  two 
gases  to  mix  before  they  reach  the  exit  plane.  The  mixture 
is  ignited  oy  arcing  a  spark  across  the  tungsten  wire  and 
the  1/2"  outer  tube.  Spark  energy  is  supplied  by  a  neon 
light  high  voltage  coil.  The  primary  coil  voltage  is  sup¬ 
plied  by  a  120  vac  outlet  followed  by  a  VARIAC  reostat. 
Spark  energy  is  controlled  by  the  VARIAC  setting  (typically 
:70vac>.  The  torch  ignitor  was  extinguished  as  soon  as  the 
flame  was  stabilized  in  the  dump  section. 

The  test  section  was  extruded  from  optical  quality 
quartz  (Supersil  )  by  Heraeus-Amer s i 1 ,  Inc.  The  cylindri¬ 
cal  test  section  has  a  152.4  mm  (6  in.)  inside  diameter 
with  3  nan  (0.120  in.)  walls  and  is  609  mm  (24  in.)  in 
length.  Although  the  quartz  itself  was  optical  quality, 
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Figure  6.  Torch  Ignitor 


I i.e.  no  bubbles,  inclusions  or  striae)  the  surface  finish 
was  not.  Circumferential  "waves"  due  to  the  extruding  pro¬ 
cess  were  present.  These  waves  had  a  period  m  the  axial 
direction  of  approximately  6.3  mm  (0-26  in.)  and  had  peaks 
and  valleys  of  0.1  mm  (:  .004  in.)  that  were  visible  to  the 
eye.  These  surface  irregularities  did  not  deviate  the  laser 
beams  to  any  major  extent  for  the  majority  of  data  points. 
However ,  it  was  necessary  to  examine  the  beam  intersection 
at  the  pinhole  of  the  PMT  to  assure  that  the  beams  were  not 
distorted.  If  they  were,  the  measurement  location  was 
shifted  slightly  to  give  a  good  quality  signal.  The  quartz 
section  was  mounted  between  two  steel  flange  plates  which 
were  compressed  by  a  set  of  four  springs.  This  allowed  for 
relative  thermal  expansion  of  the  steel  rig  and  quartz  tube 
without  the  generation  of  large  stresses.  Asbestos  rope  was 
used  as  the  gasket  material  between  the  steel  flanges  and 
the  quartz  tube.  Figure  7  shows  the  test  section  assembly. 

A  set  of  cromel-alumel  thermocouples  (Omega  #  CAIN- 
14#-12)  located  676  mm  (26.6  in.)  downstream  of  the  sudden 
expansion  allowed  the  symmetry  of  the  flow  to  be  monitored. 
The  thermocouple  beads  were  located  90°  apart  circumferen¬ 
tially.  It  should  be  noted  that  the  temperatures  measured 
were  equilibrium  temperatures  and  not  true  gas  temperatures. 
No  effort  was  made  to  shield  the  thermocouples  from  radia¬ 
tion  losses  and  no  temperature  correction  for  radiation, 
convection  or  conduction  heat  transfer  effects  was  made.  A 
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596.9 


thermocouple  rake  with  ten  thermocouples  was  used  to  measure 
the  temperature  profile  at  a  location  676  mm  downstream  of 
the  sudden  expansion.  A  ten  channel  digital  thermometer 
(Omega,  Model  2176A)  accurate  to  ±0.01°  C  was  used  to  con¬ 
vert  the  electrical  output  of  the  thermocouples  into  a  digi¬ 
tal  temperature  reading. 

The  extension  duct  was  fabricated  from  152.4  mm  (6  in.) 
standard  schedule  40  steel  pipe  and  was  approximately  1.83  m 
(6  ft.)  long.  The  length  of  this  extension  duct  was  found 
to  govern  the  range  of  fuel-air  mixture  ratios  over  which  a 
flame  could  be  stabili2ed  within  the  dump.  If  the  extension 
duct,  was  too  long  the  flame  would  be  very  unstable  with 
explosions  and  blowout  occur ing .  Shortening  the  duct 
widened  the  stability  loop  and  produced  quiet,  stable 
combust  ion 

4 .  THE  FUEL  SYSTEM 

Four  100  lb.  bottles  of  liquid  propane  connected  in 
parallel  were  used  to  supply  fuel  to  the  test  rig.  The  pro¬ 
pane  was  kept  at  33°  c  ( 90°F )  with  a  small  electrical 
heater.  This  increased  the  propane  tank  pressure  to  approx¬ 
imately  170  psig,  which  allowed  an  experimental  run  time  of 
approximately  one  hour  before  tank  pressure  dropped  signifi¬ 
cantly  [  <,  40  psig),  causing  the  run  condition  to  drift. 
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The  propane  flowed  out  of  the  tank  manifold  into  a 
regulator  set  at  40  ps  lg .  This  maintained  a  constant  flow 
rate  even  though  tank  pressure  was  decreasing  due  to  the 
latent  heat  of  vaporization.  This  regulated  flow  was  next 
split  into  two  paths.  One  path  went  to  the  torch  ignitor 
and  the  other  went  to  the  main  fuel  orifices. 

A  gas  stove  regulator  was  used  to  regulate  the  torch 
ignitor  propane  flowrate  while  a  single  stage  regulator  was 
used  to  control  the  air  flowrate.  An  electrically  con¬ 
trolled  solenoid  valve  (normally  closed)  was  placed  in  the 
propane  line  directly  upstream  of  the  torch  ignitor  for 
safety  reasons.  The  main  fuel  branch  consisted  of  a  shut¬ 
off  valve,  a  rot.ometer  (Brooks,  Model  1100),  a  two-stage 
regulator  and  an  electrically  controlled  solenoid  valve 
(normally  closed)  directly  upstream  of  the  main  fuel  mani¬ 
fold.  Figure  8  schematically  shows  the  fuel  system. 

* 

5.  THE  DATA  COLLECTION,  STORAGE  AND  PROCESSING  SYSTEM 

The  LDV  photomultiplier  tube  output  was  connected  to  a 
TS I  Model  1980  signal  processor  [49).  This  unit  features  a 
250  MHz  clock  with  two  nano-second  resolution  and  capability 
for  either  digital  or  analog  output.  Only  the  digital  out¬ 
put  was  used  in  this  investigation.  In  the  absence  of 
external  control  the  processor  data  rate  (number  of  vali¬ 
dated  velocity  measurements  per  second)  depends  on  the  rate 
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at  which  seeding  particles  enter  the  probe  volume,  the  pro¬ 
cessor  gain  setting  (which  effectively  sets  the  trigger 
level)  and  the  number  of  cycles  per  Doppler  burst  required 
for  validation.  Data  rates  as  low  as  a  few  per  second  to  in 
excess  of  20,000  per  second  were  obtained.  The  rate  at 
which  particles  enter  the  test  section  was  controlled  by  the 
particle  seeder. 

Data  acquisition  and  short  term  storage  were  performed 
by  an  IMSAI  8080  microcomputer  and  Micropolis  floppy  disk 
system.  Data  could  be  sampled  (by  microcomputer  control) 
from  rates  as  low  as  0.1  sample  per  second  to  approximately 
4800  samples  per  second.  The  microcomputer  had  the  capabil¬ 
ity  of  storing  15,600  data  points  in  its  memory. 

After  sampling,  the  data  were  written  onto  a  floppy 
disk  for  temporary  storage.  Each  disk  is  capable  of  storing 
100,000  data  points.  The  microcomputer  also  interfaced  with 
Purdue  University's  CDC  6500  and  6600  computers.  Data  may 
be  transferred  from  the  floppy  disk  to  magnetic  tape  for 
permanent  storage.  The  CDC  r-00  and  6600  were  used  for  all 
data  reduction. 

The  microcomputer  and  TS I  processor  were  also  inter¬ 
faced  so  that  velocity  sampling  was  jointly  controlled  by 
the  processor  and  the  microcomputer.  When  the  TS I  processor 
has  a  data  point  teady,  it  transmits  a  data  ready  pulse  to 
the  microcomputer.  Upon  reception  of  the  data  ready  pulse 
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cne  microcomputer  returns  a  data  inhibit  pulse  to  the  TSI 
processor.  The  inhibit  stops  the  processor  from  accepting 
more  data  and  causes  it  to  hold  the  present  data  until  it 
can  be  read.  The  microcomputer  waits  a  fixed  amount  of 
time,  chosen  by  the  operator,  before  reading  the  data  point. 
After  the  data  is  read  the  data  inhibit  is  removed  by  the 
microcomputer  and  the  cycle  continues  until  the  desired 
number  of  samples  has  been  taken.  The  rate  at  which  data  is 
actually  acquired  is  therefore  controlled  by  the  slower  of 
the  two  instruments.  When  seeding  density  is  high  and  the 
sampling  time  interval  is  large,  the  data  will  be  acquired 
at  essentially  equal  time  intervals  and  an  unbiased  velocity 
distribution  will  result  [4,5]. 

6.  THE  SEEDING  SYSTEMS 

Two  flow  seeders  were  used  in  this  experiment.  One  for 
the  isothermal  flow  and  one  for  the  reacting  flow.  The 
seeder  used  for  the  unbiased  cold  flow  was  one  which  is  com¬ 
mercially  built,  by  TSI.  It  consisted  of  a  Model  3074  air 
supply,  a  Model  3076  liquid  atomizer  and  a  Model  3072 
evaporation-condensation  monodisperse  aerosol  generator. 
This  system  produced  seeding  particles  about  one  micron  or 
less  in  diameter  using  a  solution  of  100%  Dioctyl  Phthalate 
(DOP).  The  seeder  was  operated  at  constant  pressure,  nor¬ 
mally  60  ps i ,  while  the  evaporation-condensation  unit  was 
operated  at  a  constant  voltage,  60  volts.  The  seed  density 
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inside  the  flow  section  was  sufficient  to  give  particle 
arrival  rates  in  excess  of  20,000  per  second  at  the  LDV 
probe  volume  over  most  of  the  test  region. 

Hot  flow  measurements  were  made  with  1  urn  Al^O.^  seeding 
particles  supplied  by  a  TSI  Model  3400  fluidized  bed  parti¬ 
cle  generator.  The  maximum  particle  number  density  gen¬ 
erated  by  the  aluminum  oxide  fluidized  bed  seeder  was  much 
lower  than  the  particle  density  generated  by  the  DOP  seeder. 
Particle  arrival  rates  of  5,000  per  second  in  the  central 
part  of  the  test  region  were  the  maximum  obtainable  with  the 
present  fluidized  bed  particle  generator.  The  aluminum 
oxide  seeder  was  also  used  for  making  biased  measurements  in 
the  cold  flow  so  that  a  comparison  between  cold  flow  and 
reacting  flow  results  using  the  same  type  of  seeding  could 
be  made. 
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SECTION  IV 


EXPERIMENTAL  TECHNIQUE 


1.  INTRODUCTION 

Experimental  mapping  of  the  axisymmetric  sudden  expan¬ 
sion  flow  field  consisted  of  direct  measurements  in  both 
isothermal  and  reacting  flows.  Both  biased  and  unbiased 
velocity  measurements  were  made  in  the  isothermal  flow  case, 
while  only  biased  measurements  were  made  in  the  reacting 
flow  case.  More  will  be  said  about  this  later.  This  sec¬ 
tion  presents  the  techniques  used  to  obtain  the  following 
flow  field  parameters: 


1.  mean  streamwise  velocity,  u 

2.  streamwise  turbulence  intensity,  \ I u ' ^ 

3.  stream  function,  f 

4.  reattachment  length, 

5.  integrated  mass  flow  rate 

Also  included  in  this  section  are  the  appropriate 
values  of  the  various  LDV  system  and  flow  system  parameters 
used  in  the  study. 
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2 .  TEST  PROCEDURE 


As  noted  previously,  three  sets  of  data  were  taken  at 
each  measurement  location  consisting  of  two  data  sets  for 
cold  flow  (biased  and  unbiased)  and  one  data  set  for  hot 
flow  (biased).  These  data  are  presented  in  Section  V.  Corn- 
par  isons  between  the  unbiased  and  biased  cold  flow  results 
illustrate  the  effect  of  velocity  bias  while  comparisons 
between  the  biased  hot  and  cold  flow  data  show  the  effect  of 
combustion  on  the  flow  field.  (Unbiased  hot  flow  data  could 
not  be  obtained  directly  because  of  seeder  limitations.) 

All  flow  conditions  were  maintained  at  near  constant 
values  throughout  the  testing  procedure.  The  inlet  center - 
line  velocity  was  22.07  m/s  ±.08m/s.  The  fuel  flow  was 
maintained  at  0.0038  lbm/s  ±.0002  lbm/s.  This  gives  an 
overall  fuel-air  ratio  of  0.018  and  an  overall  equivalence 
ratio  of  0.28.  The  hot  flow  run  condition  was  controlled  by 
monitoring  the  inlet  centerline  velocity  and  the  fuel  flow 
reading  on  the  rotometer.  Four  exit  thermocouples  were  also 
used  to  monitor  the  hot  flow  run  condition.  Typically,  the 
exit  temperatures  were  maintained  constant  to  within  ±10°C. 

3.  MEAN  VELOCITIES  AND  TURBULENCE  PARAMETER  MEASUREMENTS 

The  mean  velocities  and  turbulence  parameters  are  cal¬ 
culated  from  LDV  measurements  at  various  grid  points  in  the 
flow  field.  System  parameters  were  set  and  held  constant 
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throughout  the  series  of  measurements  to  insure  continuity 
of  technique  and  minimize  sources  of  error.  Table  4  is  a 
listing  of  the  system  parameters  used  for  all  measurements. 

The  techniques  employed  to  obtain  the  beam  waist  diame¬ 
ter  and  position  can  be  found  in  Ref.  [49).  The  beam  half 
angle  was  measured  and  found  to  be  3.52  ±.03  degrees. 

The  fringe  spacing,  FR,  can  be  determined  from 


F 


R 


X 

o 


1  ,  8  .  ' 

2  sin  (j) 


(1) 


0 

Substitution  of  the  measured  half  angle  ^  into  Equation 
(1)  yields  a  fringe  spacing  of  4.19  ±.04  Mm.  The  seed  par¬ 
ticle  size  was  approximately  Imih  in  diameter  for  both  the 
DOP  and  aluminum  oxide  particles. 

The  net  frequency  shift  employed,  fg,  was  10  MHz  ±1 
KHz.  The  positive  frequency  shift  indicates  that  the 
fringes  are  moving  upstream  against  the  mean  flow  direction, 
in  this  case  with  a  velocity  of  41.9  m/s.  This  was  adequate 
for  the  identification  of  negative  velocities  at  all  points 
in  the  flow  and  eliminated  incomplete  signal  bias  as  noted 
ear lier . 
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TABLE  4. 


LiH  S V S T E  M  I’AKAHLT U<  SETTINbb 


Optical  System  Settings 

XQ  (laser  frequency):  0.5145  pm 

0/2  (beam  intersection  half  angle):  3.52  °  ±  0.03  ° 

FR  (fringe  spacing):  4.19  pm  ±  .04  pm 
Probe  volume  size  (waist  diameter):  130  pm  +  5  pm 
Fs  (frequency  shift):  10  MHZ  ±  1  KH2 

Data  Collection  Settings 

Electronic  filter  settings:  30  MH  (low  pass) 

3  MH*  (high  pass) 

N  (number  of  fringes/signal):  16 

Comparator:  3  (3.1%) 

n  (exponent) :  floating 

Sample  size:  4500  samples 

Data  rate:  >  20,000  samples/sec  (unbiased) 

500  -  1500  samples/sec  (biased) 

Sample  rate:  50  samples/sec  (unbiased) 

free  running  processor  (biased) 

Seed  particles:  dioctyl  phthalate  -  D0P  (unbiased) 
aluminum  oxide  -  A^O^  (biased) 

Flow  System  Parameters 

U,  (inlet  centerline  velocity):  22.07  m/s  ±  0.08  m/s 

4 

ReH  (based  on  centerline  velocity):  5.51  x  10 

m,  (air  mass  flow  rate):  0.21  Ibm/s  t  0.01  lbm/s 

tiip  (propane  massflow  rate):  0.0038  lbm/s  ±  .0002  lbm/s 

$  (overall  equivalence  ratio):  0.28  ±  0.02 

H  (step  height) :  33. 1  mm 

R1  (inlet  radius):  33.1  mm 

R2  (outlet  radius):  76.2  mm 


The  output  signal  of  the  photomultiplier  tube  was  fil¬ 
tered  to  remove  the  "pedestal*  and  any  high  frequency  noise. 
A  30  MHz  low  pass  filter  and  a  3  MHz  high  pass  filter  were 
used.  A  16  to  8  fringe  comparison  (N  «  16)  was  used  with 
the  accuracy  of  this  comparison  set  at  3.1  percent  (compara¬ 
tor  =  3).  This  means  that  the  time  for  a  particle  to  cross 
16  fringes  is  compared  with  twice  the  time  for  the  same  par¬ 
ticle  to  cross  8  fringes.  An  error  in  the  comparison  of 
more  than  3.1  percent  results  in  the  measurement  being 
rejected.  The  exponent,  n,  was  allowed  to  float  (variable) 
and  was  an  output  of  the  TS I  processor. 

For  unbiased  cold  flow  measurements  the  liquid  atomizer 
followed  by  the  evaporation-condensation  unit  supplied  DOP 
particles  at  a  number  density  sufficient  to  give  particle 
arrival  rates  (number  of  valid  data  points  per  second  avail¬ 
able  for  sampling  by  the  TS I  processor)  in  excess  of  20,000 
per  second  at  the  probe  volume  over  most  of  the  test  region. 
Hot  flow  measurements  were  made  with  Al^O^  seeding  particles 
supplied  by  the  fluidized  bed  seeder.  Seeding  number  den¬ 
sity  was  lower  in  this  case,  with  maximum  particle  arrival 
rates  of  5,000  per  second  in  the  central  core  of  the  test 
region.  The  particle  arrival  rate  was  actually  maintained 
at  between  500  1500  particles  per  second  throughout  the 
entire  test  Beet,  ion  during  hot  flow  runs  by  varying  the  gain 
on  the  TSI  processor  and  the  laser  power.  This  was  done  to 
reduce  variability  in  the  measurements  due  to  seeding 
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density  variations.  The  fluidized  bed  seeder  was  also  used 
for  making  biased  measurements  in  the  cold  flow  case  for 
later  comparison  to  the  reacting  flow  results. 

The  method  of  velocity  bias  elimination  suggested  by 
Roesler,  et  al.  [4,5]  was  employed  in  the  cold  flow  studies 
due  to  the  high  seeding  density  obtainable  with  the  liquid 
aerosol  seeder.  Unfortunately,  the  fluidized  bed  seeder 
used  for  the  hot  flow  study  did  not  have  a  capacity  large 
enough  to  produce  the  seed  density  required  for  full  velo¬ 
city  bias  elimination.  The  processor  was  free  running  (lim¬ 
ited  only  by  the  maximum  cycling  rate  of  the  microcomputer 
in  this  case)  so  that  validated  Doppler  signals  were 
recorded  as  fast  as  they  could  be  read  into  the  microcom¬ 
puter.  The  rate  at  which  the  microcomputer  was  able  to  sam¬ 
ple  the  available  data  was  at  its  maximum  of  4800  samples 
per  second.  In  all  cases  each  velocity  measurement  was  con¬ 
structed  from  4500  samples.  The  sample  size  was  large 
enough  to  give  valid  results  without  using  excessively  large 
computer  storage.  A  maximum  validation  rate  of  1500  parti¬ 
cles  per  second  (biased  data  only)  allows  the  sampling  to 
take  place  over  a  three  second  interval.  This  interval  was 
believed  long  enough  to  contain  the  ma]or  part  of  the  tur¬ 
bulence  spectrum  in  the  present  flow  field.  Using  statisti¬ 
cal  analysis  and  assuming  gausian  distributions  and  local 
turbulence  intensity  values  of  70%  (defines  standard  devia¬ 
tion  of  histogram)  giveB  the  expected  sampling  error  for 
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this  sample  size.  Mean  velocities  and  turbulence  intensi¬ 
ties  were  found  to  be  statistically  accurate  to  ±1%  and  ±3%, 
respectively  for  a  95i  confidence  level  [50]. 

Each  sampled  output  of  the  TS I  processor  consisted  of 
three  digital  numbers  namely  N  (cycles/burst),  n  (exponent) 
and  D  (diqital  matissa).  This  data  was  converted  into  a 
frequency  f  and  then  a  velocity,  u  by  the  following  equa- 
tions  [48 ] : 
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N*  x  108 

D  x  2n-2 
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Where  f  is  the  frequency  shift  and  F  is  the  fringe  spacing 
calculated  from  Equation  (1).  The  mean,  variance  and  skew¬ 
ness  coefficient  of  the  4500  individual  velocities  were  then 
computed  from  Equations  (4),  (5)  and  (6)  below. 
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In  computing  values  from  Equations  (4),  (5)  and  (6)  any 

individual  measurements  deviating  more  than  ±3o  from  the 
mean  were  discarded  as  noise.  The  number  of  discarded 
points  was  typically  about  10  to  15  and  always  less  than  50 
per  data  set. 

The  experimental  grid  consisted  of  70  grid  points 
divided  into  7  radial  grid  lines  as  shown  in  Figure  9.  Each 
grid  point  was  spaced  approximately  7.02  mm  (0.300  in)  apart 
in  the  radial  direction  while  the  grid  line  separation  was 
two  step  heights  starting  at  x/H  =  1.  The  grid  line  at  x/H 
13  was  excluded.  The  “inlet"  velocity  and  turbulence 
intensity  profiles  were  measured  at  x/H  =  1/3.  This  meas¬ 
urement  plane  was  as  close  to  the  sudden  expansion  plane  as 
the  LDV  system  could  be  used  due  to  physical  obstruction. 

Although  the  system  geometry  was  axially  symmetric  the 
flow  downstream  of  the  sudden  expansion  may  not  be.  Thus  it 
was  necessary  to  determine  if  symmetry  existed  in  the  down¬ 
stream  flow.  Symmetry  was  tested  for  the  cold  flow  by  two 
means.  As  a  coarse  check  static  pressure  measurements  were 
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taken  circumferentially  at  six  streamwise  locations  in  a 
steel  test  section  of  the  same  geometry  as  that  of  the 
quartz  test  section.  These  measurements  showed  no  asymmetry 
in  the  flow.  The  second  more  sensitive  technique  involved 
taking  LDV  measurements  on  both  sides  of  the  centerline  of 
the  test  section.  The  results  showed  that  the  mean  flow  and 
rms  velocities  were  indeed  axially  symmetric.  Because  of 
this,  only  measurements  along  a  radius  in  the  horizontal 
plane  were  required  to  map  the  flow  field.  The  hot  flow  was 
also  believed  to  be  symmetric  as  the  fuel  was  injected  well 
upstream  of  the  sudden  expansion  at  four  circumferential 
locations.  Also,  the  amount  of  fuel  added  to  the  flow  was 
very  low  (f/a  =  0.018)  and  did  not  change  the  flow  i xeld 
appreciably. 

4.  THE  STREAM  FUNCTION 


The  stream  function,  f,  for  an  ax l symmet r 1 c  incompres¬ 
sible  flow  is  defined  by 


u 


a* 
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Using  second  order  central  differencing  about  the  points  i 
and  i  1,  Equation  (7)  can  be  aproximated  by 
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At  any  given  x-plane  F(x)  will  be  a  constant  and  can  be 
arbitrarily  set  to  zero.  Because  of  this,  Equation  (8) 
becomes 


l  +  l 


*i  + 


(r 


l  +  l 


-  r 


i) 


u  ,  .  +  u 
l  +  l  l 


(9) 


The  stream  function  values,  calculated  from  equation  (9), 
were  then  tabulated  and  contours  of  the  stream  function  were 
determined  for  specific  values  of  this  parameter. 

5.  THE  REATTACHMENT  LENGTH 

The  reattachment  length,  x^ ,  was  determined  by  linearly 
extrapolating  a  u  =  0  curve  t.o  the  wall.  The  u  *  0  curve 
was  located  by  linear  interpolation  between  adjacent  grid 
points  at  which  u  changed  sign. 


6.  THE  MASS PLOW  RATE 

Since  the  flow  was  ax lsymmetr ic ,  it  was  possible  to 
compute  the  mass  flux  at  each  measurement  plane.  This 
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permitted  a  further  check  on  the  accuracy  of  the  velocity 
data.  Both  the  biased  and  unbiased  cold  flow  mean  velocity 
data  were  integrated  using  piecewise  integration  of  a  poly¬ 
nomial  fit.  The  mass  flux  at  each  grid  line  was  normalized 
with  the  inlet  mass  flux,  which  was  calculated  using  a  1/7 
power  law  velocity  profile  for  fully  developed  pipe  flow. 
This  power  law  profile  gave  an  integrated  mass  flux  that 
agreed  with  the  experimental  integrated  mass  flux  at  grid 
lines  located  at  x/H  =  5  and  7  to  better  than  one  percent 
(see  Table  5).  Mass  flux  was  not  computed  for  the  hot  flow 
case  due  to  the  lack  of  temperature  and  therefore  density 
information. 
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SECTION  V 


EXPERIMENTAL  RESULTS 


1.  INTRODUCTION 

In  this  section  the  results  of  the  LDV  velocity  meas¬ 
urements  in  both  the  isothermal  (cold)  flow  and  reacting 
(hot)  flow  cases  are  presented.  Cold  flow  data  in  the 
biased  and  unbiased  modes  will  be  compared.  As  previously 
noted,  only  biased  hot  flow  measurements  could  be  made.  The 

Reynolds  number  for  the  axisymmetric  sudden  expansion  based 

4 

on  step  height  and  inlet  centerline  velocity  was  5.5  x  10  , 
which  corresponds  to  an  inlet  centerline  velocity, 
U^=22 . 07m/s .  Measurements  with  combustion  occur ing  were 
made  using  gaseous  propane  and  air  at  an  overall  equivalence 
ratio  of  *=0-28.  A  lean  mixture  was  employed  to  Keep  the 
temperature  in  the  test  section  at  levels  which  allowed 
steady  state  operation.  The  propane  was  injected  well 
upstream  of  the  sudden  expansion  as  noted  in  Section  III. 
Thus  the  propane  and  air  were  reasonably  well  mixed  at  the 
sudden  expansion.  However,  a  completely  premixed  mixture 
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was  not  obtained,  since  stable  combustion  was  achieved  at  an 
overall  fuel-air  mixture  ratio  well  below  the  lean  burning 
1 lmit . 

Representative  plots  of  mean  streamwise  velocity  , 
streamwise  turbulence  intensity  normalized  with  the  inlet  > 

centerline  velocity  and  the  local  streamwise  turbulence 
intensity  for  both  cold  flow  (biased  and  unbiased)  and  hot 
L low  (biased)  are  presented.  Skewness  coefficient  profiles 
at  four  grid  lines  are  also  presented.  Flow  streamlines  and 
reattachment  lengths  are  derived  from  the  experimental  data. 

Integrated  mass  flux  comparisons  of  the  biased  and  unbiased 
cold  flow  velocity  data  are  made  at  different  grid  lines  and 
the  measured  temperature  profile  at  x/H  =  17.7  is  presented. 

The  comparison  of  the  cold  flow  measurements  with  numerical 
predictions  is  contained  in  Section  VI.  The  reduced  experi¬ 
mental  velocity  data  is  tabulated  in  Appendix  Al.  The  tabu¬ 
lated  data  include  mean  velocities,  normalized  turbulence 
intensities  and  local  turbulence  intensities  for  each  meas¬ 
urement  point..  Temperature  data  at  the  exit  plane  is  tabu¬ 
lated  m  Appendix  A2 . 

2.  AVERAGE  STREAMWISE  VELOCITY 

The  inlet  velocity  profile  was  measured  experimentally 
at  x/H  =  0.33.  This  measurement  plane  was  as  close  to  the 
plane  of  the  sudden  expansion  as  the  LDV  setup  would  permit. 
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Figure  10  shows  the  experimental  unbiased  streamwise  velo¬ 
city  measurements  at  x/H  *  0.33  along  with  a  1/7  power  law 
velocity  profile.  As  the  figure  shows,  the  agreement 
between  the  measurements  and  a  1/7  power  law  profile  was 
very  good.  This  allowed  one  to  obtain  the  inlet  mass  flow 
rate  by  analytically  integrating  the  power  law  velocity  pro¬ 
file  across  the  inlet  pipe  radius.  The  power  law  profile 
was  also  used  as  an  inlet  boundary  condition  in  the  numeri¬ 
cal  prediction  code  described  in  Section  VI. 

Figures  11  through  17  show  the  normalized  mean  Btream- 
wise  velocity  profiles  at  measurement  planes  located  1,  3, 
5,  7,  9,  11,  and  15  step  heights  downstream  of  the  sudden 
expansion.  The  data  sets  presented  in  each  of  the  figures 
are  for  biased  and  unbiased  cold  flow  and  for  biased  hot 
flow  Figure  18  summarizes  the  normalized  mean  strearawise 
velocity  data  presented  in  Figures  11  through  17  (cold  flow 
only)  m  one  plot.  The  computational  predictions  described 
in  Section  VII  are  also  shown  on  this  figure.  Figure  19 
summarizes  biased  hot  and  unbiased  cold  flow  normalized  mean 
streamwise  velocity  data  in  one  plot.  Although  the  statist¬ 
ical  error  was  found  to  be  rather  small  (±1%  for  mean  velo¬ 
cities  and  ±3%  for  turbulence  intensities),  there  can  be 
additional  experimental  error  due  to  the  problem  of  defining 
the  probe  volume  location  exactly  and  run  condition  drift. 
The  experimental  error  can  be  identified  roughly  by  noting 
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the  deviation  between  individual  data  points  and  the  solid 
line  drawn  through  the  data.  The  solid  lines  represent  a 
"best  fit"  profile. 

Several  observations  can  be  made  from  the  mean  stream- 
wise  velocity  data.  The  trends  of  the  biased  versus 
unbiased  mean  velocities  are  as  one  would  expect.  That  is, 
the  absolute  value  of  the  biased  mean  velocity  is  higher 
than  the  absolute  value  of  the  unbiased  mean  velocity 
throughout  the  flow  field.  The  differences  are  largest  in 
the  regions  of  high  turbulence  intensity  and  are  insignifi¬ 
cant  in  regions  of  low  turbulence  intensity.  There  is  a 
crossover  of  the  mean  velocity  profiles  as  you  pass  from  a 
positive  flow  region  into  a  negative  flow  region  as  shown  in 
Figures  11  through  13.  Typically,  the  biased  cold  flow  data 
was  less  than  1  m/s  higher  than  the  unbiased  cold  flow  data 
throughout  the  flow  field.  Velocity  differences  in  the 
biased  and  unbiased  cold  flow  data  are  almost  indistinguish¬ 
able  in  the  core  of  the  inlet  jet,  but  show  a  difference  of 
from  0.4  to  0.6  m/s  in  the  turbulent  region  at  axial  dis¬ 
tances  greater  than  five  step  heights. 

Maximum  negative  velocities  in  the  recirculation  zone 
were  -2.9  m/s  and  -4.8  m/s  for  the  cold  flow  and  hot  flow 
cases,  respectively.  These  values  convert  to  -13%  and  -22% 
of  the  inlet  centerline  velocity.  Although  the  recircula¬ 
tion  zone  was  found  to  be  somewhat  stronger  in  the  hot  flow 
case  it  was  also  found  to  be  thinner  and  shorter  than  the 
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cold  flow  recirculat ion  zone. 


Pitz  also  observed  this 


effect  in  Ref.  (29]  where  LDV  measurements  were  made  in 
isothermal  and  reacting  flows  behind  a  rearward  facing  step. 
The  profiles  in  Pigure  19  show  that  reattachment  occurs 
between  7  and  9  step  heights  downstream  of  the  sudden  expan¬ 
sion  in  the  cold  flow  case  and  at  approximately  7  step 
heights  downstream  of  the  sudden  expansion  in  the  hot  flow 
case.  The  effects  of  combustion  (i.e.  higher  velocities  due 
to  higher  temperatures)  are  not  seen  until  axial  distances 
are  greater  than  three  step  heights.  Note  the  inflection 
point  present  in  the  velocity  profiles  at  x/H  =  15  in  Pigure 
17.  This  character ist ic  was  also  found  by  Stevenson,  et  al. 
(15)  and  indicates  that  the  flow  is  far  from  fully 
developed. 

Pigure  20  shows  the  centerline  velocity  decay  for  the 
unbiased  cold  flow  and  the  biased  hot  flow  cases.  Also 
shown  in  this  figure  are  the  corresponding  cold  flow  data 
obtained  by  Stevenson,  et  al.  (15],  Moon  and  Rudinger  (12) 
and  Preeman  (11J.  The  cold  flow  data  obtained  in  the 
present  study  show  very  good  agreement  with  that  of  Steven¬ 
son,  et  al.  (15],  Moon  and  Rudinger  (12)  and  Freeman  (11) 
when  plotted  as  a  function  of  downstream  distance  normalized 
with  step  height.  This  is  in  spite  of  the  fact  that  the 
expansion  area  ratios  of  the  various  test  sections  were  dif¬ 
ferent  by  as  much  as  a  factor  of  two.  Also,  inlet  boundary 
conditions  were  different,  varying  from  an  almost  flat  inlet 
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Figure  20.  Measured  Mean  Centerline  Velocity  Decay 
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velocity  profile  in  (15],  to  fully  developed  pipe  profiles. 
This  figure  also  indicates  that  the  flow  has  not  relaxed  to 
fully  developed  pipe  flow  at  IS  step  heights  downstream  of 
the  sudden  expansion,  since  CJ/U1  *  0.41  instead  of  0.25  as 
would  be  obtained  from  an  area  ratio  argument. 


Measurement  inaccuracies  due  to  steep  velocity  gra¬ 
dients  can  be  a  problem  when  using  an  LDV  system.  The 
errors  in  high  velocity  gradient  regions  can  be  large,  espe¬ 
cially  if  there  are  fluctuations  (instabilities)  in  the  flow 
field.  This  error  is  primarily  due  to  the  finite  probe 
volume  size.  For  the  geometry  used  in  this  study,  the 
highest  velocity  gradients  occur  in  the  shear  layer  at  a 
radius  equal  to  the  inlet  pipe  radius  close  to  the  sudden 
expansion  plane.  Using  the  velocity  profile  at  x/H  »  1  to 
calculate  the  velocity  gradient  across  the  shear  layer  gives 
6u/6R«140s  This  converts  to  a  velocity  gradient  of  0.14 
m/s  pet  millimeter.  In  this  study  the  probe  volume  (length 
of  2  mm)  would  therefore  have  a  0.28  m/s  velocity  gradient 
across  it.  This  does  not  effect  the  mean  velocity  calcula¬ 
tion,  but  does  effect  the  turbulence  intensity  calculation. 
Using  the  procedure  given  by  Karpuk  and  Tiederman  (51), 
which  assumes  a  cylindrical  probe  volume  and  a  linear  velo- 
city  gradient  across  the  probe  volume,  the  maximum  error  in 
turbulence  intensity  due  to  the  finite  size  of  the  probe 
volume  was  found  to  be  0.08  m/s.  Thus  velocity  gradient 
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error  was  deemed  insignificant  for  the  present  investiga¬ 
tion,  since  this  calculation  is  a  worst  case  example. 


3.  STREAMWISE  TURBULENCE  INTENSITY 

The  normalized  inlet  turbulence  intensity  profile  at 
x/H  =  0.33  is  shown  in  Figure  21.  As  expected  fairly  low 
levels  of  turbulence  are  present  in  the  central  core  of  the 
inlet  "jet"  flow,  whereas  large  levels  of  turbulence  are 
present  in  the  shear  layer  at  the  edge  of  the  inlet  pipe. 
This  figure  shows  that  the  turbulence  intensity  normalized 
to  the  inlet  centerline  velocity  varies  from  3.5%  in  the 
core  region  to  15%  in  the  shear  layer.  Since  the  local 
average  velocity  decreases  as  you  move  toward  the  shear 
layer  the  local  inlet  turbulence  intensity  varies  from  3.5% 
at  the  centerline  to  31.3%  in  the  shear  layer  at  x/H  =  0.33. 

Pigures  22  through  28  show  the  streamwise  turbulence 
intensity  normalized  with  the  inlet  centerline  velocity  at 
measurement  planes  located  1,  3,  5,  7,  9,  11,  and  15  step 
heights  downstream  of  the  sudden  expansion  plane.  Again, 
results  from  the  three  data  sets  mentioned  earlier  are 
presented  here.  Pigure  29  sunmarizes  the  normalized  stream- 
wise  turbulence  intensity  data  presented  in  Pigures  22 
through  28.  These  figures  all  show  a  peak  in  the  normalized 
turbulence  intensity  which  broadens  as  the  mixing  zone 
spreads  in  the  downstream  direction.  The  peak  turbulence 
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intensities  identity  regions  ot  nign  snear  arid  mixing.  Max¬ 
imum  normalized  streamwise  turbulence  intensities  were  found 
to  be  approximately  22%  for  both  the  hot  and  cold  flows.  By 
comparison,  maximum  values  of  normalized  turbulence  inten¬ 
sity  reported  for  backward  facing  step  flows  are  28%  by  Pitz 
{29],  21%  by  Eaton  and  Johnston  [5,6]  and  19%  by  Brenmer,  et 
al.  (3J.  These  results  are  in  good  agreement  with  the 
present  study.  Downstream  of  the  recirculation  zone,  at  x/H 
=■  15,  the  st teamwise  turbulence  intensity  profile  was  found 
to  be  almost  flat  with  normalized  turbulence  intensities  of 
approximately  15%. 

One  interesting  result  found  in  the  cold  flow  study  was 
the  change  in  shape  of  the  normalized  turbulence  intensity 
profiles  as  a  result  of  velocity  bias  as  shown  in  Pigures  22 
through  27.  These  changes  in  profile  were  found  to  be 
insignificant  in  the  low  turbulence  core  region  where  bias¬ 
ing  effects  were  small,  and  significant  in  regions  of  high 
turbulence.  Most  notably,  the  location  of  maximum  shear  was 
found  to  differ  significantly  for  the  biased  and  unbiased 
data  sets.  The  value  of  normalized  turbulence  intensity 
obtained  at  a  given  location  in  a  highly  turbulent  region 
was  also  found  to  differ  (by  as  much  as  35%).  It  should  be 
noted  that  these  differences  were  due  only  to  the  sampling 
technigue,  as  the  fluid  structure  of  the  flow  remained  the 
same.  This  indicates  that  the  probability  distributions 
constructed  from  the  individual  velocity  realizations  were 
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controlled  not  only  by  the  turbulence  in  the  flow,  but  by 
the  sampling  technique.  These  trends  agree  with  the  trends 
of  the  biased  and  two-dimensional  corrected  data  of  Tieder- 
man  (47]  as  previously  shown  in  Pigure  3.  This  indicates 
that  the  two-dimensional  correction  scheme  appears  to  change 
the  velocity  probability  distribution  toward  an  unbiased 
one. 

Because  biased  data  gives  a  higher  mean  velocity  than 
unbiased  data  (more  high  velocity  particles  pass  through  the 
probe  volume  per  unit  time  than  low  velocity  particles),  one 
would  expect  biased  data  to  give  lower  local  turbulence 
intensities  than  unbiased  data.  This  is  because  the  vari¬ 
ance,  as  defined  by  Equation  5,  is  divided  by  a  larger  local 
mean  velocity  in  the  biased  case  as  opposed  to  the  unbiased 
case.  This  assumes,  of  course,  that  there  is  no  change  in 
the  variance  due  to  the  sampling  technique.  In  the  outer 
regions  of  the  flow  the  biased  mean  velocities  were  found  to 
be  approximately  1  m/s  higher  (in  positive  flow  regions) 
than  the  unbiased  mean  velocities  as  mentioned  in  Section 
4.2.  Although  the  variance  of  the  biased  measurements  was 
found  to  be  larger  (higher  normalized  turbulence  intensity) 
in  this  same  region,  the  effect  of  the  higher  local  mean 
velocities  (biased)  was  to  reduce  the  local  turbulence 
intensity  for  biased  measurements.  This  overriding  effect 
was  because  the  1  m/s  mean  velocity  offset  due  to  bias  was  a 
ma^ot  portion  of  the  absolute  local  mean  velocity  in  this 


outer  region.  Figures  30  through  36  show  the  local  stream- 
wise  turbulence  intensity  profiles  for  the  biased  and 
unbiased  cold  flow  on  a  semi-log  scale  and  verifies  this 
result.  It  should  be  noted  that  values  of  local  turbulence 
intensity  greater  than  100-200%  can  be  deceiving,  because 
they  reflect  regions  of  near  zero  mean  velocity.  Typically, 
m  non-zero  mean  velocity  regions,  the  biased  local  tur¬ 
bulence  intensity  is  25-35%  less  than  the  unbiased  local 
turbulence  intensity.  In  the  potential  core  region,  at 
axial  distances  less  than  five  step  heights,  the  difference 
was  found  to  be  insignificant  as  expected. 

It  is  interesting  to  note  that  the  centerline  local 
turbulence  intensity  increases  rapidly  as  you  move  down¬ 
stream  from  the  potential  core.  Figure  37  shows  the  local 
centerline  turbulence  intensity  for  the  unbiased  cold  flow 
data  from  this  study  along  with  data  from  Stevenson,  et  al. 
(15].  Good  agreement  was  obtained  considering  the  different 
inlet  boundary  conditions  used  in  the  two  studies.  The 
local  turbulence  intensity  increases  from  14%  for  0  <  x/H  < 
5  to  a  maximum  value  of  38%  at  x/H  =  15. 

This  again  indicates  the  complex  turbulent  structure  in 
the  region  downstream  of  the  potential  core.  Local  values 
of  centerline  turbulence  intensity  in  a  fully  developed  pipe 
flow  are  approximately  4%.  The  values  obtained  here,  how¬ 
ever,  are  of  the  same  magnitude  as  those  found  in  the  near 
wall  region  of  turbulent  pipe  flow. 
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Figure  30.  Measured  Local  Streamwise  Turbulence  Intensity  Profiles  at  x/H 
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Figure  31.  Measured  Local  Streanwise  Turbulence  Intensity  Profiles  at  x/H 
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Figure  33.  Measured  Local  Streamwise  Turbulence  Intensity  Profiles  at  x/H 
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igure  34,  Measures  Local  Streamwise  Turbulence  Intens 
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Figure  35.  Measured  Local  Streamwise  Turbulence  Intensity  Profiles  at  x/H 


Comparisons  between  the  biased  hot  and  cold  flow  nor¬ 
malized  turbulence  intensities  show  that  the  recirculation 


zone  is  narrower  and  shorter  in  the  hot  flow  case.  This 
result  was  also  found  in  the  two-dimensional  rearward  facing 
step  study  by  Pitz  [29].  Maximum  normalized  turbulence 
intensities  for  the  hot  flow  case  were  approximately  22%. 
This  is  the  same  maximum  value  found  in  the  cold  flow  case. 
It  is  interesting  to  note  that  the  maximum  normalized  tur¬ 
bulence  intensity  at  x/H  =  1  (Figure  22)  for  the  hot  flow 
case  was  33%  higher  than  the  value  found  for  the  cold  flow 
case.  Pitz  [29]  also  found  that  the  reacting  flow  had  nor¬ 
malized  turbulence  intensities  30-35%  higher  than  the  isoth¬ 
ermal  flow  at  this  location.  The  reason  for  this  is  not 
c lear . 

Although  the  normalized  turbulence  intensities  were 
found  to  be  of  the  same  magnitude  for  cold  and  hot  flow,  the 
local  turbulence  intensities  were  not.  Figures  30  through 
36  show  that  the  hot  flow  local  turbulence  intensities  are 
much  lower,  typically  40  to  60%  of  the  cold  flow  values,  at 
planes  located  greater  than  three  step  heights  downstream  of 
the  sudden  expansion.  Pitz's  data  [29]  also  exhibited  this 
result.  Combustion  has  not  affected  the  flow  field  at  planes 
upstream  of  three  step  heights,  especially  in  the  central 
core  teg  ion.  At  the  plane  located  15  step  heights  down¬ 
stream  of  the  sudden  expansion  the  local  values  of  stream- 
wise  turbulence  intensity  were  found  to  vary  from  38%  on  the 


center  line  to  774  at  the  neat  wall  point  in  the  cold  flow 
case,  whereas  the  corresponding  values  m  the  hot  flow  case 
were  74  and  514.  Pigure  37  also  shows  that  hot  flow  local 
centerline  turbulence  intensity  was  much  less  than  in  the 
cold  flow. 

The  fact  that  local  turbulence  levels  in  the  reacting 
flow  were  lower  than  in  the  isothermal  flow  may  be  because 
dilatation  by  heat  release  competed  with  or  dominated  tur¬ 
bulence  energy  production  by  shear  flow,  notably  in  the 
vicinity  of  the  recirculation  boundary.  This  result  was 
also  found  by  Pu")ii  and  Eguchi  |22)  in  their  study  where  UDV 
measurements  were  made  in  a  bluff  body  flame  stabilizer. 


4.  SKEWNESS  COEFFICIENTS 

The  skewness  coefficient  is  a  statistical  parameter 
commonly  used  to  help  identify  the  turbulence  structure.  It 
is  a  measure  of  the  symmetry  or  asymmetry  of  the  velocity 
probability  distribution  and  is  calculated  by  using  Equation 
6.  By  definition  biased  velocity  distributions  are  more 
negatively  skewed  (a  greater  number  of  high  velocity  reali¬ 
zations)  than  unbiased  velocity  distributions  .  This  was 
indeed  found  to  be  the  case  as  shown  in  the  skewness  coeffi¬ 
cient  profiles  plotted  in  Pigures  38  and  39.  Note  that  for 
values  of  r/R2  less  than  0.55,  the  skewness  coefficient 
(which  indicates  distribution  shape)  is  similar  for  both  the 
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unbiased  and  biased  measurements .  This  would  indicate  that 
the  shapes  of  the  histograms  constructed  from  the  two  sam¬ 
pling  techniques  were  similar.  That  is,  the  effect  of  bias¬ 
ing  in  the  central  core,  where  low  turbulence  exists,  is 
minimal.  However,  as  one  moves  radially  outward  toward  the 
point  where  the  normalized  turbulence  intensity  profiles 
cross,  the  distributions  start  to  deviate  significantly  from 
one  another. 

The  skewness  of  the  reacting  flow  measurements  looks 
quite  different  from  that  in  the  cold  flow  at  first  glance. 
However  there  are  similarities  when  one  considers  that  the 
reacting  shear  layer  is  shifted  closer  to  the  test  section 
wall.  This  is  the  case  before  reattachment  occurs  and  can 
be  seen  at  x/H  =  3  and  7  in  Figures  38  and  39,  respectively. 
If  the  shift  in  the  shear  layer  location  for  hot  flow  versus 
cold  flow  is  taken  into  account  Figures  38  and  39  show  that 
shapes  of  the  histograms  are  very  similar,  although  non¬ 
gauss  lan,  at  x/H  »  3  and  7.  This  was  not  the  case  at  planes 
downstream  from  reattachment.  as  shown  in  Figures  40  and  41. 
The  skewness  was  much  more  negative  in  the  reacting  flow 
case  indicating  different  turbulence  structure.  The  skew¬ 
ness  tends  to  be  highly  negative  at  the  centerline  and  is 
zero  or  slightly  positive  at  the  wall  in  the  reacting  flow. 
Since  no  unbiased  hot  flow  data  were  taken,  it  is  not  known 
how  much  of  this  difference  in  structure  (skewness  coeffi¬ 
cient)  was  due  to  the  combustion  process  and  how  much  was 


due  to  sampling  technique,  although  the  small  difference  in 
skewness  observed  between  the  biased  and  unbiased  cold  flow 
data  would  seem  to  indicate  that  combustion  was  the  dominant 
factor . 

5.  THE  REATTACHMENT  LENGTH 

The  reattachment  point  was  determined  by  linearly 
extrapolating  the  u*0  contour  to  the  wall  boundary.  The 
measured  reattachment  length  for  the  unbiased  cold  flow  was 
found  to  be  8.62  step  heights  downstream  of  the  step  face. 
This  value  is  in  good  agreement  with  earlier  results 
reported  in  the  literature  111,12,13,15).  The  measured 
reattachment  length  for  the  biased  hot  flow  was  found  to  be 
7.40  step  heights  downstream  of  the  step  face.  The  reat¬ 
tachment  length  for  the  reacting  shear  layer  was  reduced  by 
'.15%.  Pitz  found  reductions  in  reattachment  lengths  between 
20  and  30%  in  his  study. 


6.  THE  STREAM  PONCTION 

Figures  42  and  43  show  lines  of  constant  volume  flow 
for  the  unbiased  cold  flow  and  the  biased  hot  flow,  respec¬ 
tively.  These  are  stream  function  contours  for  the  cold 
flow  (constant  density)  case.  To  partially  compensate  for 
slight  mass  flow  discrepancies  due  to  experimental  errors, 
the  stream  function  was  non-dimensionalized  at  each  x/H 
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station.  That  is  at  each  x  station 


* 


* 

*wall . 


(10) 


The  unbiased  cold  flow  measurements  showed  that  the 
center  of  the  recirculation  region  was  located  3.5  step 
heights  downstream  of  the  step  face  at  a  radial  distance  of 
r /R 2  =  0.725.  The  center  of  the  recirculating  region  in  the 
hot  flow  case  was  located  at  x/H  *  3.0  and  r/R2  *  0.750. 
This  again  illustrates  the  general  shift  of  the  shear  layer 
toward  the  wall  in  the  hot  flow  case. 


7.  THE  MAS SF LOW  RATE 

In  this  study  two  sets  of  streamwise  velocity  measure¬ 
ments  for  the  isothermal  flow  were  taken.  As  pointed  out 
earlier  the  first  set  of  measurements,  considered  isothermal 
flow  were  taken.  The  first  set  of  measurements,  considered 
to  be  unbiased,  was  taken  with  particle  arrival  rates  in 
excess  of  20,000  Hz  and  a  computer  sampling  rate  of  50  sam¬ 
ples  per  second.  The  second  set  of  data  was  taken  with  a 
particle  arrival  rate  between  500  and  1500  Hz  and  free  sam¬ 
pling,  that  is,  the  microcomputer  was  allowed  to  sample  data 
as  fast  as  it  could  (4800  Hz).  These  two  mean  velocity  data 
sets  were  then  reduced  and  integrated  using  piecewise 
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integration  with  forward 
points  and  backward  direct 
f  it. 


polynomial  fitting  to  the  data 
integration  of  the  polynomial 


Table  5  lists  the  normalized  integrated  mass  flux  at 
each  axial  measurement  plane  for  the  unbiased  and  biased 
measurements.  The  mass  flux  was  normalized  with  the  inlet 
mass  flux  which  was  obtained  by  analytically  integrating  a 
1/7  power  law  velocity  profile.  The  measured  values  indi¬ 
cate  that  the  unbiased  integrated  mass  fluxes  agree  to 
within  ±3%  from  plane  to  plane.  The  biased  measurements 
overpr edict  normalized  integrated  mass  flux  by  as  much  as 
24%.  This  result  again  demonstrates  the  effect  of  velocity 
bias  and  gives  one  confidence  in  the  unbiased  velocity  meas¬ 
urement  technigue  used.  Unfortunately,  integration  of  mass 
flux  in  the  reacting  flow  was  not  possible  as  detailed  tem¬ 
perature  (density)  information  was  not  available. 


8.  TEMPERATURE  PROFILES 

A  temperature  traverse  in  the  reacting  flow  was  made  at 
a  plane  located  17.7  step  heights  downstream  of  the  sudden 
expansion.  Figure  44  shows  the  measured  stagnation  tempera¬ 
ture  contours.  From  the  contours,  one  can  see  that  the  core 
region  was  fairly  cool  which  indicates  that  the  mixing  pro¬ 
cess  was  far  from  complete  at  this  location.  The  tempera¬ 
tures  here  were,  as  pointed  out  earlier,  not  true  gas 
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TABLE  5. 

INTEGRATED  MASS  1  LUX 
SUDDEN  EXPANSION 

IN  AN  AX  1  SYMMETRIC 

x/H 

Biased  Data 

Unbiased  Data 

iii/m 

re  f 

■n/mref 

0.33 

1.07 

1.00* 

1.00 

1.07 

1.04 

3.0 

1.20 

1.06 

5.0 

1.13 

1.00 

7.0 

1.21 

1.01 

9.0 

1.24 

1.03 

11.0 

1.21 

1.02 

15.0 

1.12 

1.03 

*mref  based  on  biased  inlet  velocity  profile  at  x/H  =  0.33 


1 20 


Top 

(Looking  Downstream) 


Figure  44.  Temperature  Contours  at  x/H 
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temperatures  but  "equilibrium"  temperatures,  because  no 
effort  was  made  to  correct  the  thermocouple  readings  for 
heat  transfer  losses.  The  main  purpose  for  this  temperature 
traverse  was  to  determine  profile  symmetry.  Pigure  44  shows 
that  the  temperature  profile  was  synsneti  ic  across  the  meas¬ 
urement.  axis  but  was  not  symmetric  from  top  to  bottom.  The 
reason  for  this  asymmetry  is  not  known,  but  it  may  be  due  to 
fuel  infection  imbalances.  An  attempt  was  made  to  improve 
symmetry  by  adjusting  the  flow  in  the  four  fuel  inlets,  but 
Figure  44  represents  the  best  result  obtainable.  The  asym¬ 
metry  was  believed  not  severe  enough  to  cause  significant 
errors  in  the  velocity  data  taken  on  the  horizontal  central 
plane,  however. 
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SECTION  VI 

COMPARISON  OP  NUMERICAL  ANALYSIS 
WITH  ISOTHERMAL  EXPERIMENTAL  RESULTS 


1.  INTRODUCTION 

In  order  to  provide  a  basis  for  compat  ison  of  analyti¬ 
cally  predicted  and  experimentally  measured  isothermal  flow 
parameters,  the  computer  code  CHAMPION  2/E/P IX  of  Pun  and 
Spalding  [52]  was  adapted  to  the  flow  geometry  and  run. 
This  code  uses  the  k  -  t  turbulence  model  of  Launder  and 
Spalding  [IS]  along  with  a  modified  version  of  the  SIMPLE 
algorithm  described  by  Pantankar  and  Spalding  [52].  Unfor¬ 
tunately,  a  prediction  code  capable  of  modeling  the  combus¬ 
tion  process  was  not  available  for  use  in  the  present  study. 
Therefore  only  isothermal  comparisons  will  be  made  in  this 
sect  ion . 

The  2/E/PIX  code,  when  adapted  for  an  axtsymmetr ic 
geometry,  solves  the  partial  differential  equation 


a 

3x 


(pu«fr)+i 


d 

tFr 


(prvc<J»)  - 


„  d 

Va* 


I'.Sl 


,  i 

r 


d 

37 


dd> 

(vf?> 


(io; 


123 


where  r  is  the  coordinate  in  the  radial  direction,  x  is  the 


coordinate  in  the  atreamwise  direction,  u  is  the  mean  velo¬ 
city  in  the  atreamwise  direction,  vf  1a  the  mean  velocity  in 
the  radial  direction,  <t>  is  the  dependent  variable,  is  the 
exchange  coefficient,  and  s^  is  the  source  term.  Equation 
(10)  represents  the  time  averaged  Navier-Stokes  equations 
along  with  a  supplemental  transport  equation.  The  k  -  e 
model  assumes  isotropic  diffusion  with  the  effective  viscos¬ 
ity,  being  the  sum  of  the  laminar  and  turbulent  con¬ 
tributions.  That  ia. 


eff 


^lam  + 


t. 


(11) 


When  appropriate  expressions  for  4>,  sAand  r  as  listed  in 
Table  6  are  substituted  into  the  general  Equation  (10),  the 
equation  takes  on  the  form  of  continuity,  axial  and  radial 
momentum,  turbulent  kinetic  energy  and  rate  of  energy  dissi¬ 
pation  equations. 

The  2/E/FIX  code  solves  this  set  of  s imultaneous  dif¬ 
ferential  equations,  with  appropriate  boundary  conditions, 
using  upwind  differencing.  The  matrix  equation  obtained 
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TABLE  6.  CONSERVATION  EQUATIONS  CORRESPONDING  TO  EQUATION  10 


from  the  numerical  approximations  is  then  solved  using  a 
tridiagonal  algorithm  along  with  under  relaxation  to  achieve 
numerical  stability. 

Inputs  to  the  code  include  inlet  velocity,  inlet  tur¬ 
bulent  kinetic  energy,  relaxation  factors,  and  five  con¬ 
stants  used  in  the  expressions  listed  in  Table  6.  The 
numerical  values  of  these  five  constants  recommended  by 
Launder  and  Spalding  [20]  are  given  m  Table  7. 

TABLE  7.  RECOMMENDED  TURBULENCE  CONSTANTS  FROM  REP.  [20] 


Although  the 
Table  6,  it 


istant 

Value 

CD 

0.09 

C1 

1.43 

C2 

1.92 

°K 

1.00 

1.30 

is  in 

none  of 

is  used  in  the  following  two 


the  expressions 
equations 


in 


^turb 


c 


(12) 
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Equation  12  is  the  equation  from  which  the  turbulent 
contribution  to  viscosity  is  evaluated.  Equation  13 
represents  a  modified  "log  law"  used  to  link,  the  first  node 
from  the  wall  to  the  wall  boundary  condition  where  k  is  a 
constant  equal  to  0.4,  E  is  a  constant  equal  to  9.0  for 

smooth  walls,  r  is  the  wall  shear  stress,  u  is  the  mean 

w 

axial  velocity  at  point  y^,  and  y^  is  the  distance  from  the 
wall  to  the  first  node  from  the  wall. 

A  21  x  41  grid  was  used  fot  the  computations.  The  21 
radial  increments  covered  the  tube  radius  while  the  41  axial 
grid  planes  were  located  at  each  integer  value  of  x/H  (down¬ 
stream  distance  normalized  to  step  height)  from  x/H  =  0  to 
40.  At  x/H  =  40  the  flow  is  essentially  fully  developed  and 
an  exit  boundary  condition  of  zero  axial  velocity  gradient 
is  valid. 

The  velocity  profile  measured  as  close  as  physically 
possible  to  the  plane  of  the  sudden  expansion  (x/H  =  0.33) 
was  used  as  the  inlet  boundary  condition.  The  convergence 
criterion,  built  into  the  code,  on  u,  v ,  turbulent  kinetic 
energy  (TKE)  and  turbulent  energy  dissipation  (TED)  was 
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2.  MATCHING  RE ATTACHMENT  LENGTH 

Reattachment  length  was  chosen  as  a  common  flow  field 
parameter  for  comparison  of  the  numerical  analysis  to  exper¬ 
imental  results-  Moon  and  Rudinger  [12]  matched  their 
experimental  results  in  terms  of  reattachment  length  by  set¬ 
ting  the  coefficients  C^=1.43  and  C2=l-70.  In  the  study  by 
Stevenson,  et  al .  [15],  only  was  iterated  upon  to  match 
the  experimental  value  of  reattachment  length.  A  value  of 
C2=l-94  was  used  in  that  study.  In  the  present  study  C 2  was 
also  the  only  coefficient  iterated  upon.  A  least  square 
linear  fit  was  then  applied  to  the  data,  yielding  the  fol¬ 
lowing  equation. 


x 

C2  -  0.059  +  2.4305  (14) 

where  x^  is  the  reattachment  length  and  H  is  the  step 
height.  Since  the  observed  reattachment  length  was  approxi¬ 
mately  8.6  step  heights,  a  value  of  1.92  was  chosen  for 
based  on  this  equation.  This  is,  in  fact,  the  value  sug¬ 
gested  by  Pun  and  Spalding. 


3.  MEAN  STREAMWISE  VELOCITY 

Figures  45  through  51  show  streamwise  mean  velocity 
profile  comparisons  of  the  experimental  unbiased  cold  flow 


128 


Figure  45.  Comparison  of  Predicted  and  Measured  Mean  Streamwise  Velocity  Profile  at  x/H 
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Figure  47.  Comparison  of  Predicted  and  Measured  Mean  Streamwi 
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Figure  49.  Comparison  of  Predicted  and  Measured  Mean  Streamwise  Velocity  Profile  at  x/H 
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Figure  51.  Comparison  of  Predicted  and  Measured  Mean  Streamwi 


data  and  the  numerically  predicted  results  at  different 
axial  locations.  Figure  52  shows  the  normalized  stream 
function  cor.cours.  It  can  be  seen  that  the  computational 
results  are  in  fairly  good  agreement  with  the  unbiased  data, 
although  some  discrepancies  exist  in  the  recirculation  zone 
and  the  predicted  velocities  in  the  center  of  the  flow  are 
consistently  low  at  downstream  locations.  This  agrees  with 
earlier  results  [15].  The  figures  also  show  that  as  the 
flow  develops  the  profiles  begin  to  deviate  more.  Although 
there  is  a  large  difference  in  the  experimental  and 
predicted  velocities  in  the  center  of  the  flow,  the  differ¬ 
ence  in  mass  flow  rate  is  small  since  cross-sectional  flow 
area  in  this  region  is  small.  tn  order  to  check  the  code, 
the  integrated  mass  flux  was  calculated  using  the  predicted 
velocities  at  each  grid  line.  Mass  conservation  was  satis¬ 
fied  experimentally  by  the  unbiased  measurements  as  shown 
earlier  in  Table  5,  but  the  computer  code  did  not  conserve 
mass.  Typically,  3.21  of  the  inlet  mass  flux  was  "lost"  by 
the  x/H  =  15  grid  line.  This  was  true  for  all  values  of 
tested.  The  mass  flux  loss  was  attributed  to  truncation  and 
convergence  errors  in  the  finite  difference  scheme.  The 
method  used  to  link  the  pressure  field  with  the  velocity 
field  may  have  also  caused  slight  errors. 
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SECTION  VII 


CONCLUSIONS  AND  RECOMMENDATIONS 

Several  conclusions  can  be  drawn  from  the  results  of 
this  investigation.  Measurements  of  the  mean  streamwise 
velocity  field  in  the  presence  of  combustion  showed  the 
anticipated  changes  relative  to  the  cold  flow  case.  Those 
changes  include  highet  mean  velocities  due  to  heat  release 
and  a  shorter,  stronger  recirculation  zone.  The  effects  of 
combustion  were  not  influential  at  planes  upstream  of  three 
step  heights. 

Turbulence  intensity  measurements  indicated  significant 
differences  in  shear  layer  position  and  rate  of  turbulence 
decay  for  the  hot  flow  case,  but  relatively  small  differ¬ 
ences  in  peak  normalized  turbulence  levels.  Maximum  normal¬ 
ized  turbulence  intensities  were  found  to  be  approximately 
22%  in  both  the  hot  and  cold  flows,  while  local  turbulence 
intensities  were  found  to  be  much  lower  in  the  hot  flow 
case,  typically,  40  to  60%  of  the  cold  flow  values.  This 
seems  to  indicate  that  the  turbulence  production  mechanism 
due  to  shear  is  similar  in  both  flow  cases,  but  that  combus¬ 
tion  tends  to  suppress  the  turbulence  in  regions  away  from 
the  shear  layer . 
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The  results  obtained  in  this  study  have  again  verified 
the  fact  that  the  velocity  bias  elimination  method  based  on 
high  seeding  density  and  fixed  interval  data  sampling  is 
effective  for  highly  turbulent  internal  flows.  Normalized 
mass  flux  values  calculated  using  unbiased  mean  velocity 
data  agreed  throughout  the  flow  field  to  within  ±3%,  while 
the  calculations  using  the  biased  mean  velocity  data  over- 
predicted  the  mass  flux  by  as  much  as  24*.  The  sampling 
technique  itself  was  also  found  to  "bias"  the  statistical 
turbulent  quantities  (i.e.  turbulence  intensity  and  skew¬ 
ness)  siqnif icantly  in  highly  turbulent  legions. 

Numerical  prediction  of  the  mean  velocity  field  using 
the  CHAMPION  2/E/FIX  code  was  reasonably  accurate,  although 
computed  velocities  in  the  center  of  the  flow  were  substan¬ 
tially  below  those  measured  at  downstream  locations.  A 
ma^or  deficiency  in  this  code  was  the  lack  of  mass  flow  con¬ 
servation  which  may  have  led  to  the  poor  prediction  of  mean 
centerline  velocities. 

Some  areas  which  need  further  investigation  are  as  fol¬ 
lows  : 

1)  unbiased  reacting  flow  measurements  to  determine 
directly  the  velocity  bias 

2)  an  effective  velocity  bias  correction  or 
elimination  scheme  that  could  be  applied 

to  flow  mesurements  in  sparsely  seeded  flow 
fields  including  those  with  significant 
density  fluctuations 

3)  LDV  measurements  of  tangential  and  radial  mean 
velocities  and  turbulence  intensities  along  with 
Reynolds  stress  correlations  in  both  isothermal 
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and  reacting  flows;  knowing  these  values  and 
using  the  Favre  averaged  turbulent  kinetic 
energy  equation  for  this  type  of  flow  could 
explain  the  combustion-turbulence  interaction 
in  more  detail 

4)  the  effect  of  axial  pressure  gradient  on  the 
turbulence  structure  of  the  flow  field  also 
seems  to  be  important  and  its  effects  need  to 
be  investigated  systematically. 

As  a  final  comment  it  should  be  noted  that  no  special 
problems  were  encountered  when  using  the  LDV  in  this  combus¬ 
tion  study.  No  observable  degradation  in  signal  quality  was 
caused  by  refractive  index  fluctuations  in  the  flow. 
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EXPERIMENTAL  DATA 
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Table  AT.  Experimental  Velocity  and  Turbulence 
Data  ( x/H* 1 ) 


R 

R/R2 

D 

u/ui 

/u'/ui 

(mm) 

(m/s) 

v 

COLO 

FLOW  (UNBIASED) 

0 

0 

21.963 

.995 

.038 

.038 

7.620 

.100 

21.701 

.933 

.044 

.045 

15.240 

.200 

20.661 

.936 

.053 

.057 

22.860 

.300 

19.332 

.876 

.063 

.072 

30.480 

.400 

17.060 

.773 

.075 

.097 

38.100 

.500 

11.048 

.501 

.  151 

.301 

45.720 

.600 

.390 

.018 

.089 

5.049 

53.340 

.700 

-.612 

-.028 

.057 

-2. OSS 

60.960 

.800 

-.616 

-.028 

.067 

-2.403 

68.580 

.900 

-.411 

-.019 

.075 

-4.010 

COLD  FLOW  (BIASED) 

0 

0 

22.015 

.998 

.038 

.038 

7.620 

.100 

21.892 

.992 

.045 

.045 

15.240 

.200 

20.777 

.941 

.054 

.058 

22.860 

.300 

19.641 

.890 

.061 

.069 

30.480 

.400 

17.361 

.787 

.076 

.096 

38.100 

.500 

11.523 

.522 

.142 

.27! 

45.720 

.600 

2.129 

.096 

.127 

1.319 

53.340 

.700 

-1.064 

-.048 

.070 

-1.452 

60.360 

.800 

-1.013 

-.046 

.067 

-1.463 

68.580 

.900 

-.021 

-.037 

.073 

-1.970 

HOT 

FLOW  (BIASED) 

0 

0 

22.254 

1.008 

.040 

.040 

7.620 

.100 

22.221 

1.007 

.042 

.041 

15.240 

.200 

21.446 

.972 

.051 

.052 

22.860 

.300 

19.886 

.901 

.062 

.069 

30.480 

.400 

18.026 

.817 

.074 

.090 

38.100 

.500 

14.E29 

.663 

.  102 

.153 

45.720 

.600 

3.810 

.173 

.146 

.844 

53.340 

.700 

.904 

.041 

.202 

4.928 

60.960 

.800 

-.924 

-.042 

.084 

-2.002 

68.580 

.900 

-.250 

-.011 

.068 

-6.001 

7/E/FIX 

0 

0 

22.575 

1.023 

- 

- 

3.810 

.050 

22.575 

1.023 

- 

- 

7.620 

.100 

22.155 

1.004 

- 

- 

11.430 

.150 

21.735 

.985 

- 

- 

15.240 

.200 

21.210 

.961 

- 

- 

19.050 

.250 

20.560 

.932 

- 

- 

22.860 

.300 

19.740 

.894 

- 

- 

26.670 

.350 

18.690 

.847 

- 

- 

30.480 

.400 

17.010 

.771 

- 

- 

34.290 

.450 

14.385 

.652 

- 

- 

38.100 

.500 

8.673 

.393 

- 

- 

41.910 

.550 

5.145 

.233 

- 

- 

45.720 

.600 

3.202 

.  145 

- 

49.530 

.650 

1.302 

.  059 

- 

- 

53.340 

.700 

-.109 

-.005 

- 

- 

57.150 

.750 

-.930 

-.042 

- 

- 

60.960 

.800 

-1.428 

-.065 

- 

- 

64.770 

.850 

-1.785 

-.081 

- 

- 

68.580 

.300 

-2.079 

-.094 

- 

72.390 

.950 

-2.415 

-.109 

- 
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Table 

A1 .  cont 

.  Experimental 

Velocity 

and 

Turbulence  Data  (x/H= 

3) 

/7£ 

ca , 

R 

R/R2 

u 

U/'Ul 

W  U  SUI 

■sju  /U 

(mm) 

Wj) 

COLO 

FLOW  (UNBIASED) 

0 

0 

22.0*10 

.359 

.041 

.041 

7.G20 

.  10  0 

22. 112 

1.002 

.044 

.044 

15.2<VJ 

.200 

20.713 

.333 

.054 

.053 

22.8S0 

.300 

13.397 

.073 

.071 

.080 

20.480 

.000 

18.90; 

.770 

.114 

.143 

33.100 

.500 

10.652 

.403 

.170 

.352 

•15.720 

.COO 

5.  149 

.83.3 

.165 

.711 

53.3*10 

.700 

.640 

.089 

.  137 

4.G62 

60.960 

.200 

-1.869 

-.050 

.  100 

-1.705 

83.530 

.300 

-2.301 

-.  131 

.003 

-.674 

COLO 

FLOW  (BIASED) 

0 

0 

22.100 

1.001 

.033 

.033 

7.620 

.100 

21.893 

.979 

.045 

.046 

15.200 

.200 

?  0.7*!  4 

.040 

.054 

.037 

22.880 

.300 

13. -107 

.833 

.  OC-5 

.074 

30.080 

.400 

16.477 

.74  7 

.108 

.144 

30.100 

.500 

1 1 .90G 

.533 

.162 

.300 

05.720 

.600 

6.771 

.307 

.153 

.519 

53.300 

.700 

3.202 

.  t43 

.149 

1.003 

60.960 

.800 

-1.203 

-.055 

.139 

-2.547 

GO. 530 

.300 

-2.335 

-.136 

.094 

-.689 

HOT  FLOW  (BIASED) 

0 

0 

22.275 

1.003 

.042 

.042 

7.620 

.100 

22.471 

1.013 

.040 

.033 

15.200 

.200 

81.597 

.373 

.053 

.054 

22.860 

.300 

20.307 

.980 

.063 

.063 

30.030 

.400 

1C. 354 

.G38 

.073 

.083 

30.100 

.500 

15.130 

.623 

.112 

.162 

05.720 

.300 

10.786 

.4(16 

.1*10 

.288 

53.300 

.700 

3.C24 

.173 

.  165 

.955 

GO. 360 

.800 

-2.605 

-.113 

.132 

-1.119 

68.530 

.300 

->1.819 

-.218 

.002 

-.376 

2/E/FlX 

0 

0 

21.040 

.390 

- 

- 

3.010 

.030 

8!  .640 

.950 

- 

- 

7.620 

.  too 

21.3-5 

.  975 

- 

- 

11.030 

.  150 

81.000 

.052 

- 

- 

15.2*10 

.200 

20.4.5 

.928 

- 

- 

13.050 

.850 

19.740 

.094 

- 

- 

22.660 

.300 

U'.C.MO 

.847 

- 

- 

c w . 670 

.383 

17.010 

.771 

- 

- 

30.080 

.400 

14.700 

.ccr. 

- 

- 

30.290 

.*150 

11.9/0 

.543 

- 

- 

30.100 

,500 

3.  l'  G 

.415 

- 

- 

01.310 

.553 

e.a-t; 

.3"  3 

- 

05.720 

.('00 

>i.rfn 

.203 

- 

- 

49.S30 

.  C50 

8.5-  6 

.130 

- 

S3. 300 

.700 

1  ••ill.' 

.06-1 

- 

- 

57.150 

.  730 

.  uu 

.  oou 

- 

- 

GO. 550 

.COO 

— .  1  >  i  6 

-.0-10 

- 

- 

6*1.770 

.650 

-1.753 

-.073 

- 

- 

66.580 

.300 

-.>.499 

-.11.1 

- 

- 

72.390 

.510 

-3. -183 

155 

- 

143 


Table  Al.cont.  Experimental  Velocity  and 
Turbulence  Data  (x/H=5) 


COLD  FLOW  (UNBIASED) 


0 

0 

21 .763 

.  00G 

.015 

.046 

7.620 

.100 

21.434 

.371 

.051 

.052 

15.240 

.200 

20.841 

.917 

.  0G7 

.073 

22.(560 

.300 

17.360 

.737 

.113 

.151 

30.480 

.400 

14.563 

.860 

.163 

.247 

32. 100 

.500 

8.066 

.452 

.130 

.420 

05.720 

.600 

4.581 

.208 

.  107 

.303 

53.340 

.700 

.014 

.001 

.150 

231.103 

G0.8SO 

.800 

-.501 

-.023 

.144 

-B.343 

GS.580 

.300 

-1.351 

-.008 

.11G 

-1.308 

COLD 

FLOW  (BIASED) 

0 

0 

21.368 

.935 

.045 

.045 

7.620 

.  100 

21.560 

.977 

.051 

.052 

15.240 

.200 

20.404 

.325 

.  0G4 

.063 

22.660 

.300 

10. 181 

.084 

.100 

.122 

30.4C0 

,400 

16. 144 

.731 

.140 

.131 

38.100 

.500 

11.474 

.520 

.172 

.330 

45.720 

.600 

7.018 

.318 

.104 

.578 

53.340 

.700 

.810 

.037 

.135 

5.252 

GO. 360 

.800 

-1.562 

-.072 

.  1G3 

-2.2G3 

GO. 580 

.300 

-2.000 

-.091 

.127 

-1.400 

HOT 

FLOW  (BIASED) 

0 

0 

22.162 

1.004 

.042 

.042 

7.620 

.100 

21.SC7 

.956 

.044 

.044 

15.240 

.200 

21.533 

.976 

.054 

.055 

<'2 . 660 

.300 

20 . 855 

.  31U 

.066 

.072 

30.480 

.400 

18.119 

.021 

.081 

.093 

r.r.  too 

.500 

1 4.675 

.  61 .5 

.113 

.130 

45.720 

.600 

10.382 

.435 

.147 

.287 

53.310 

.700 

7.005 

.317 

.168 

.533 

80.860 

.800 

1.089 

.086 

.225 

2.623 

r.n.sao 

.300 

-4.372 

-.  r.)3 

.  135 

-.680 

2/E/FIX 

0 

0 

20.665 

.937 

- 

- 

3.810 

.  050 

20.685 

.337 

- 

- 

7.620 

.100 

20.370 

.323 

- 

- 

11.430 

.150 

10.740 

.834 

- 

- 

15.240 

.200 

13;900 

•  85G 

- 

- 

10.050 

.850 

17.745 

.804 

- 

- 

22.860 

.300 

16.065 

.  7811 

- 

- 

86.670 

.350 

14.175 

•  G42 

- 

- 

30.480 

.400 

18. 180 

.352 

- 

- 

34.230 

.450 

10.090 

.457 

- 

- 

33. 100 

.300 

0.13/ 

•  3G9 

- 

- 

41.310 

.  550 

6.363 

.280 

- 

- 

45.720 

.600 

4.767 

.216 

- 

- 

43.330 

.650 

3.360 

.  152 

- 

- 

53.340 

.700 

2.110 

.096 

- 

- 

57. 150 

.760 

1 . 0  1 H 

•  046 

- 

- 

60 . 560 

.800 

.066 

.003 

- 

- 

64.770 

.8.30 

- .  760 

-.034 

- 

- 

68.580 

.8ii0 

-1 .4  >0 

-.0G7 

- 

- 

72.330 

.'.150 

-8.343 

-.  ion 

- 

- 

144 


Table 

A1 .  cont 

.  Experimental  Velocity 

and 

Turbulence  Data 

( x/H  = 

7) 

/1,2 

/T!i 

R 

R/R2 

u 

0/4)1  •> 

J  u  /ui 

a/U  /U 

(mm) 

(m/s) 

v 

COLO 

FLOW  (UNBIASED) 

0 

0 

13.915 

.902 

.  053 

.070 

7.320 

.  100 

13.345 

.904 

.075 

.033 

15.240 

.800 

18.377 

.833 

.110 

.141 

22.RR0 

.300 

15.435 

.899 

.171 

.244 

no.-ino 

.400 

11.132 

.504 

.198 

.390 

33. ICO 

.500 

8.235 

.378 

.205 

.545 

■15.720 

.800 

4. 843 

.210 

.  106 

.884 

53.340 

.700 

1.817 

.082 

.187 

2.025 

80.950 

.800 

.882 

.030 

.151 

5.035 

63.5R0 

.300 

-.711 

-.032 

.123 

-3.827 

COLO 

FLOW  (BIASED) 

0 

0 

20.851 

.345 

.054 

.087 

7.620 

.100 

20.445 

.928 

.070 

.078 

15.240 

.200 

19.755 

.895 

.031 

.102 

22.G60 

.300 

18.531 

.749 

.148 

.198 

30. 480 

.400 

12.810 

.580 

.178 

.307 

38.100 

.500 

9.233 

.418 

.  179 

.428 

45.720 

.800 

8.722 

.305 

.203 

.885 

53.3-10 

.700 

3.081 

.139 

.203 

1 .432 

RO.2S0 

.000 

1.188 

.053 

.175 

3.308 

E3.580 

.300 

-.839 

-.029 

.158 

-5.462 

HOT  FLOW  (BIASED) 

0 

0 

21.S34 

.978 

.051 

.052 

7.820 

.100 

21.748 

.985 

.049 

.050 

15.840 

.200 

20.550 

.932 

.055 

.059 

22.360 

.300 

19.997 

.908 

.085 

.071 

30.430 

.400 

18.974 

.759 

.095 

.124 

37!.  100 

.500 

14.039 

.535 

.  134 

.210 

45.720 

.800 

10.307 

.457 

.  152 

.326 

53.340 

.700 

7.455 

.338 

.187 

.494 

60.960 

.800 

3.  133 

.  142 

.210 

1.478 

en.sao 

.SOO 

-.330 

-.017 

.178 

-10.209 

2/E/FIX 

0 

0 

19.005 

.051 

- 

- 

3.810 

.050 

19.005 

.851 

- 

- 

7.820 

.100 

18.375 

.333 

- 

- 

11.430 

.  150 

17.325 

.71)5 

- 

- 

15.340 

.800 

16.170 

.733 

- 

- 

13.050 

.250 

14.805 

.371 

- 

- 

28.880 

.300 

13.230 

.539 

- 

- 

28.870 

.330 

11.835 

.530 

- 

- 

30.480 

.400 

10. 182 

.453 

- 

- 

34.850 

.480 

8.531 

.391 

- 

- 

33.100 

.500 

7.284 

.387 

- 

- 

41.310 

.350 

5.922 

.380 

- 

- 

45.720 

.  800 

4.789 

.214 

- 

- 

40.550 

.850 

3.843 

.  185 

- 

- 

53.340 

.700 

8.557 

.  121 

- 

- 

57. 130 

.754 

l .  795 

.0111 

- 

- 

3 0.330 

.800 

1 . 020 

.047 

- 

- 

84.770 

.  850 

.351 

.018 

- 

- 

P.8.580 

.300 

-.853 

-.Oil 

- 

- 

78.530 

.350 

- 1 . 029 

-.017 

- 

- 

14b 


Table  Al.cont.  Experimental  Velocity  and 
Turbulence  Data  (x/H=9) 


R 

£✓02 

0 

Q^ui 

fn>/s) 

COLO 

flow  (unbiased) 

0 

0 

17. COO 

.797 

.134 

.168 

2.68o 

.030 

17.533 

.797 

.  133 

.173 

lo.a-'.o 

.  134 

16.633 

.754 

.160 

.212 

17.080 

.234 

1  1.991 

.679 

.132 

.268 

85.400 

.334 

12.095 

.540 

.203 

.370 

29.  100 

.430 

8.753 

.397 

.205 

.517 

■10. 720 

.534 

7.037 

.319 

.202 

.633 

on.  3*10 

.630 

4.442 

.201 

.130 

.933 

S5.9G0 

.730 

2.350 

.106 

.162 

1.525 

S3 . 530 

.034 

.865 

.039 

.139 

3.552 

COLO 

FLOW  (BIASED) 

0 

0 

13.300 

.829 

.082 

.099 

2.620 

.034 

19.272 

.828 

.104 

.126 

lo.a^o 

.134 

17.G33 

.799 

.  125 

.157 

1 7 . 360 

.234 

15.202 

.689 

.  163 

.245 

25.080 

.334 

18.831 

.584 

.104 

.316 

33. 100 

.434 

10.214 

.463 

.134 

.399 

40.720 

.534 

8.  184 

.371 

.135 

.527 

00.300 

.634 

6.  185 

.280 

.102 

.650 

55.300 

.734 

3.530 

.  1P2 

.mo 

1.149 

03.530 

.334 

1.647 

.  075 

.176 

2.355 

HOT  FLOW  (BIASED) 

0 

0 

80. POO 

.942 

.041 

.043 

2.  GOO 

.034 

80 . 785 

.939 

.042 

.044 

tO. 200 

.134 

80.046 

.  903 

.053 

.059 

17.GC0 

.214 

18.831 

.85G 

.060 

.080 

25.000 

.334 

17.216 

.  780 

.  094 

.121 

33. 100 

.434 

15.550 

.706 

.115 

.163 

00.720 

.534 

18.557 

.  5t>3 

.140 

.245 

08.300 

.634 

9.210 

.417 

.  1G7 

.400 

55.330 

.734 

5.564 

.252 

.  160 

.666 

63.580 

.334 

3.7GG 

.172 

.174 

1.015 

2/E/FIX 

0 

0 

15.645 

.709 

3.810 

.050 

15.645 

.709 

_ 

7 .  F20 

.100 

15.015 

.680 

- 

11.030 

.150 

14.070 

.633 

- 

13.200 

.200 

19. 125 

.  595 

_ 

in. 050 

.250 

11.970 

.542 

_ 

22.800 

.300 

10.815 

.  490 

_ 

26.670 

.350 

9.733 

.441 

_ 

_ 

30.080 

.000 

U.G20 

.391 

- 

30.290 

.450 

7.549 

.342 

_ 

_ 

38.100 

.500 

6.531 

.296 

- 

_ 

01.310 

.580 

5.5P6 

.853 

- 

05.720 

.600 

1 . 704 

.8)3 

_ 

03.530 

.650 

3 .  cr.5 

.  176 

. 

59.300 

.  700 

3.150 

.  143 

- 

_ 

57.150 

.750 

8.1 73 

.  112 

_ 

60.360 

.800 

1  .879 

.095 

- 

60.770 

.850 

1 .351 

.  r.«:  1 

_ 

68.580 

.300 

,•863 

•  040 

- 

_ 

72.390 

.550 

.256 

.012 

- 

- 

146 


Table  Al.cont.  Experimental  Velocity  and 
Turbulence  Data  (x/H=ll) 


K 

R'R2 

0 

Q.'Ut 

/Ta 

VU  /ui 

Vd  'V 

(mm) 

(m/t) 

COLO 

FLOW  (UNBIASED) 

0 

0 

15.016 

.63  0 

.  mn 

.291 

7.020 

.100 

M.6'10 

.  G63 

.  137 

.297 

15.240 

.200 

13. *137 

.603 

.203 

.333 

22.  rso 

.300 

11.761 

.510 

.214 

.4  19 

30.400 

.000 

H.U77 

.<102 

.203 

.504 

sn. 100 

.500 

fi.325 

.314 

.  135 

.620 

<15.720 

•  •  >o0 

6.210 

.836 

.  1713 

.753 

53.340 

.  700 

3.715 

.  147 

.  156 

1.059 

cn.nGO 

.600 

2.073 

.  038 

.133 

1.446 

60.560 

.300 

1 .265 

.057 

.118 

2.0CG 

COLD 

FLOW  (BIASED) 

0 

0 

16.042 

.772 

.163 

.228 

r.aso 

.100 

15.374 

.  697 

.177 

.254 

15.240 

.200 

13.355 

.607 

.103 

.305 

22.1180 

.300 

11.0-13 

.537 

.  103 

.351 

30. -100 

.<100 

10.043 

.455 

.  132 

.421 

3n.ioo 

.500 

8.713 

.355 

.  187 

.473 

on. 720 

.GOO 

6 . 236 

.  263 

.101 

.640 

S3. 300 

.700 

4.603 

.203 

.172 

.826 

£0.960 

.600 

3.330 

.151 

.161 

1 .065 

Co. 520 

.300 

1.463 

.057 

.133 

2.071 

HOT 

FLOW  (BIASED) 

0 

0 

20.081 

.910 

.040 

.053 

7.820 

.100 

20 .  C.CO 

.336 

.04? 

.050 

15. 200 

.  200 

13. 184 

.  L’63 

.056 

.064 

22. CEO 

.300 

18.711 

.835 

.  07  6 

.032 

3o.'teo 

.<’00 

16.768 

.737 

.105 

.142 

32.100 

.500 

13.022 

.681 

.116 

.170 

‘is,  720 

.COO 

13.014 

.530 

.  125 

.213 

53.340 

.700 

10. 135 

.453 

.  14? 

.320 

60.5G0 

.800 

7.3-35 

.351 

.  153 

.440 

63.530 

,500 

5. -US 

.245 

.111 

.453 

2/e/nx 

0 

0 

1 1 . 760 

.533 

- 

- 

3. mo 

.050 

11.760 

.5  J3 

- 

- 

7.620 

.100 

1 1.445 

.513 

- 

- 

11.030 

.150 

10.320 

.<135 

- 

- 

15.200 

.700 

10.230 

.465 

- 

- 

l '1.050 

.250 

s. ciu 

.•156 

- 

- 

72.260 

.300 

6.CT3 

.403 

- 

- 

72,670 

.  330 

153 

.370 

- 

- 

30. -100 

.<100 

7.413 

.  336 

- 

- 

31.230 

.  <150 

6.673 

.303 

- 

- 

33. 1O0 

.300 

5.555 

.872 

- 

- 

<11.310 

.550 

5.383 

.241 

- 

- 

'!*».  770 

.roo 

4.C33 

.213 

- 

- 

<15.530 

.650 

4. 116 

.  Lir, 

- 

- 

53. 3<l0 

.700 

3.570 

.  182 

- 

- 

57. ISO 

.  750 

4.076 

.133 

- 

- 

60. SCO 

.600 

8.T8  j 

.113 

- 

- 

r,  1.770 

.1  30 

8.786 

.101 

- 

- 

65.660 

.  ::oo 

t.ro'i 

.  o'ri 

- 

- 

77.350 

.  350 

1.3(3 

.  068 

- 

- 

147 


Table  Al.cont.  Experimental  Velocity  and 
Turbulence  Data  (x/H=15) 


R 

(mm) 

FVR2 

0 

(m/s) 

ii-  U 1 

yr//ui 

v/uP>u 

COLD 

FLOW  IUN8IASED) 

0 

0 

3.3GR 

.  -124 

.  1G3 

.384 

7.620 

.100 

9.15G 

.-115 

.171 

.412 

ls.a-io 

.200 

3.297 

.  42 1 

.  165 

.393 

22.1:30 

.300 

8. 456 

.295 

.  162 

.421 

30..  mo 

.400 

7.502 

.  340 

.  15-1 

.454 

38. 100 

.500 

6.468 

.294 

.  1-19 

.507 

‘IS.  720 

.eoo 

•1.974 

.  225 

.  189 

.572 

53.3-10 

.700 

4 . 039 

.184 

.117 

.63G 

Ro.cr.o 

.coo 

3.353 

.  153 

.116 

.753 

68.530 

.Duo 

2.914 

.132 

.  102 

.773 

COLD 

FLOW  (BIASED) 

0 

0 

9.311 

.449 

.  160 

.35G 

7.F20 

.100 

9.945 

.451 

.  154 

.3-11 

15.2-10 

.200 

9.673 

.  4  39 

.  160 

.385 

22. CEO 

.300 

8 . 754 

.397 

.  154 

.387 

30 .  -ISO 

,  400 

7.  13G 

.323 

.  14G 

.452 

33. 100 

.500 

6.930 

.314 

.  142 

.453 

■15.720 

.GOO 

5.377 

.244 

.120 

.493 

53.3-10 

.700 

•1.926 

.223 

.113 

.507 

SO. 530 

.800 

4.440 

.201 

.  108 

.535 

nn.sso 

.900 

3.900 

.177 

.  104 

.590 

HOT 

FLOW  (BIASED) 

0 

0 

13.342 

.876 

.  0G2 

.071 

7.620 

.  100 

19.625 

.  669 

.070 

.079 

15.2-10 

.200 

15.330 

.831 

.  063 

.105 

22.660 

.300 

16.393 

.770 

.112 

.145 

30.-::’0 

.400 

14.7G3 

.IT, 9 

.  122 

.182 

33.  10:1 

.500 

12.949 

.537 

.  168 

.213 

•15.720 

.600 

11.799 

.  535 

.127 

.237 

53. 3'!  0 

.700 

10.761 

•  463 

.  120 

.246 

SO. 830 

.500 

P .  942 

.405 

.114 

.231 

nn.500 

.300 

5.173 

.235 

.  120 

.510 

2/E/FlX 

0 

0 

7.093 

.322 

- 

“ 

3. RIO 

.050 

7 . 093 

.322 

- 

7 .620 

.mo 

6.993 

.317 

- 

- 

11.-130 

.150 

6 . 65G 

.311 

- 

“ 

15.2<:0 

.200 

6i  678 

.303 

- 

- 

in. oso 

.250 

6.4C0 

.293 

- 

- 

22.FC0 

.300 

6.237 

.  783 

- 

- 

?r,.r,70 

.330 

5.985 

.871 

- 

30.-150 

.400 

5.722 

.259 

- 

3-1.250 

.  -130 

5. '139 

.246 

- 

“ 

29. IPO 

.500 

3.  ICR 

.234 

- 

- 

•11.510 

.550 

•1.652 

.221 

- 

“ 

•15.720 

.600 

4.609 

.20.) 

- 

- 

45.530 

.630 

4. 347 

.197 

- 

- 

53.3-10 

.700 

•1.053 

.116 

- 

- 

57. 150 

.750 

3.66-1 

.  1  75 

- 

- 

ro.cno 

.non 

3.633 

.  16  1 

- 

r.4.770 

.650 

3.433 

.  156 

- 

65.530 

.500 

3.234 

.*•17 

- 

“ 

72.330 

.550 

2.930 

.  Ill 

- 

- 
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Table  A2 .  Experimental  Temperature  Data 


r/R 

_ _ 

THERMOCOUPLE  RAKE  LOOKING  DOWNSTREAM* 

TOP 

LEFT 

BOTTOM 

RIGHT 

0.9 

705 

820 

834 

861 

0.8 

682 

809 

785 

868 

0.7 

630 

785 

721 

848 

0.6 

552 

671 

595 

750 

0.5 

471 

560 

482 

647 

0.4 

397 

441 

385 

527 

0.3 

322 

350 

313 

414 

0.2 

275 

293 

275 

332 

0.1 

250 

260 

246 

271 

0.0 

232 

235 

232 

235 

♦Thermocouple  rake  located  at  x/H=17.7  and  all 
temperatures  are  in  degrees  Centigrade 
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